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I.  INTRODUCTION 


This  volume  is  part  II  of  the  final  report,  "Solid  State  Array 
Evaluation  Relevant  to  OTP  Regulations,"  [1]#  The  final  report  is  con- 
cerned with  the  evaluation  of  current  microwave  solid  state  devices 
and  an  assessment  of  their  potential  as  array  elements  with  regard  to 
the  regulations  of  the  Office  of  Telecommunications  Policy.  The  regu- 
lations are  being  revised  and  are  not  yet  in  final  form;  we  refer  to 
the  current  regulations  [24]. 

In  brief,  the  questions  posed  are  "Wlich  of  the  microwave  solid 
state' devices  are  viable  candidates  for  array  elements  in  view  of  the 
OTP  regulations  on  in-band  and  out-of-band  emissions"  and,  "What  are 
the  device  requirements  which  result  from  OTP  regulations  on  emissions?" 

Part  I of  the  final  report  [1]  is  concerned  with  an  experimental 
evaluation  of  some  current  microwave  solid  state  devices  (L  and  S-band 
bipolar  transistor  modules  and  TRAPATT  devices  were  evaluated).  Out-of- 
band  characteristics  and  the  effect  of  varying  load  were  studied  extensively. 

Part  II  (the  present  volume)  of  the  final  report  is  concerned  with 
the  computation  of  random  effects  in  planar  arrays.  Random  effects  in 
arrays,  in  particular  random  errors  in  phase  and  amplitude  excitation  of 
the  elements  of  an  array,  can  significantly  affect  the  beam  pattern  charac- 
teristics of  an  array  and  thus  affect  whether  or  not  the  OTP  regulations 
are  satisfied. 


* 


These  random  effects  are  not  extremely  significant  at  fundamental 
frequency,  for  well-designed  modules  at  L band,  for  instance.  Typical 
variations  from  module  to  module  are  + 1 db  in  amplitude  and  +10?  in 
phase  under  matched  conditions.  However,  at  harmonic  frequencies, 
significant  variations  in  amplitude  and  phase  are  encountered,  as  noted 
in  part  I [1].  In  addition,  one  obtains  larger  variations  at  the  funda- 
mental under  conditions  of  variable  load. 

The  variations  noted  in  part  I [1]  are  significant  and  must  be 
taken  into  account  in  array  analysis.  Effects  of  random  phase  and  ampli- 
tude variations  have  been  studied  in  the  past  [2-7].  These  earlier  studies 
do  not  take  mutuals  into  account  (It  turns  out  that  the  effect  of  mutuals 
is  significant.). 

A computer  program  has  been  prepared  for  the  analysis  of  planar 
arrays  of  thin-wire  dipoles  with  random  excitation  errors.  The  method 
of  moments  [8,9]  is  used  in  conjunction  with  an  iterative  method  and 
mutuals  are  taken  fully  into  account.  Figure  1 shows  the  physical  model 
studied,  a planar  array  of  dipoles  with  uniform  spacing  in  the  x,z  direc- 
tions. The  array  is  driven  by  voltage  sources  with  a nominal  (specified) 
amplitude  and  phase  distribution.  Random  amplitude  and  phase  errors  may 
be  specified.  Average  beam  patterns  are  computed  and  their  range  of 
variation  is  indicated.  The  array  may  be  treated  at  harmonic  Cor  other) 
frequencies  and  thus  an  average  ratio  of  harmonic  to  fundamental  levels 
may  be  obtained  ( as  a function  of  space  coordinates).  Expected  gain 
levels  may  also  be  computed.  A number  of  steps  have  been  taken  to  improve 
computational  efficiency  and  thereby  permit  the  treatment  of  planar  arrays 
of  moderate  size.  The  block-Toeplitz  redundancies  of  the  impedance  matrix 


and  Che  zeros  of  the  excitation  matrix  are  utilized  in  a special,  effi- 
cient solution  routine.  Arrays  of  about  100  dipole  elements  may  be 
treated  at  frequencies  up  to  about  the  third  harmonic  (i.e.,  at  fre- 
quencies where  each  dipole  is  of  length  3A/2).  This  involves  the  solu- 
tion of  a system  of  700  unknowns. 

The  results  of  computation  show  that,  with  large  phase  and  ampli- 
tude errors,  narrow  beams  are  formed  at  varying  positions  (0  ,0  ) in 

o o 

space,  with  gain  reduced  significantly  from  that  of  a uniformly  excited 
array  (with  no  errors).  For  nominal  broadside  excitation  and  completely 
random  phase,  it  is  shown  that  there  is  a tendency  to  focus  the  beams 
near  broadside,  in  contrast  to  analysis  results  without  mutuals.  The 
approximate  variation  of  gain  quantities  and  array  size  is  indicated. 

The  statistical  formulation  is  presented  in  chapter  two.  The  theory 
on  which  the  computation  is  based  is  outlined  in  chapter  three,  and  a 
computer  program  description  is  given  in  chapter  four.  Some  typical 
results  are  presented  in  chapter  five.  Tae  application  of  the  OTP  regula- 
tions is  discussed  in  chapter  six,  and  conclusions  are  presented  in 
chapter  seven.  Some  details  of  the  matrix  inversion  routine  are  given 
in  Appendix  A and  a computer  program  listing  is  given  in  Appendix  B, 

Our  conclusions  are  in  brief  that  well-designed  bipolar  transistor 
modules  and  TRAPATT  devices  are  reasonable  candidates  for  array  elements 
(assuming  that  reasonable  array  design  procedures  are  used  to  ensure 
required  side  lobe  levels) , 


4 


CL tS&SB 


2.  STATISTICAL  FORMULATION 

2.1  THE  GENERAL  STATISTICAL  FORMULATION  FOR  THE  METHOD  OF  MOMENTS 

In  the  analysis  of  random  effects  in  antenna  arrays , there  are 
essentially  many  factors  to  be  considered  as  sources  of  randomness, 
for  example,  positions  of  the  array  elements,  characteristics  of  the 
array  elements,  amplitude  and  phase  variations  of  the  excitation  vol- 
tages, etc.  Many  efforts  have  been  made  previously  to  find  the  effects 
of  these  random  errors  upon  far  field  beam  patterns  or  gain  [2-7]. 

In  this  report,  it  is  assumed  that  the  randomness  of  the  array 
is  restricted  solely  to  the  random  errors  present  in  the  voltage  exci- 
tations. The  question  is:  given  the  statistics  of  the  inputs  (voltage 

excitations)  what  will  be  the  statistics  of  the  outputs  (far  fields)? 

In  order  to  solve  the  linear  equations  of  an  antenna  array,  a 


procedure  called  the  matrix  method,  or  method  of  moments,  has  been 
used.  The  basic  mathematical  concept  is  the  method  of  moments  [8] 

by  which  the  functional  equations  of  field  theory  are  reduced  to  matrix 
equations. 

Consider  deterministic  equations  of  the  inhomogeneous  type 

F (e)  = v (2-1) 

where  F is  a linear  operator  over  linear  functional  space  V,  v is  known, 
e is  to  be  determined.  Let  e be  expanded  in  a series  of  functions 


e. , e. , e ....  in  the  domain  of  F,  as 
1 2 e 


where  a are  constants.  The  e are  called  expansion  functions.  For 
n n 

the  exact  solution,  Eq.(?-1)  is  usually  an  infinite  summation  and  e^ 
forms  a complete  set  of  basis  functions.  Substituting  Eq.  (2-2)  into 
Eq.  (2-1)  and  using  the  linearity  of  F,  one  has 


l a F(e  ) = v 
n n 
n 


(2-3) 


It  is  assumed  that  a suitable  inner  product  <e,v>  has  been  determined. 

Define  a set  of  weighting  functions,  w^,  w0 , w^,  ...  in  the  range  of  F, 

and  take  the  inner  product  of  Eq.  (2-3)  with  each  w . The  result  is 

m 


7 a <w  ,Fe  > = <w  , v> 
**  n m n m 

n 


(2-4) 


m = 1,  2,  3, 


where 


This  set  of  equations  can  be  written  in  matrix  form  as 


[f  „][<*]  = [VJ 

mn  n n 


<W1 »Fel>  <W1 tFe2> . . . 

<w2  »Fei>  <w2  »Fe2> 


(2-5) 


(2-6) 


[a  ] = a 

n i. 


[v  ] = <w  ,v> 
Q I 1 

<w2 ,v> 


(2-7) 


If  the  matrix  [f]  is  non-singular,  its  inverse  [f  ] exists.  The  cx^  are 


then  given  by 


= iCnHvl 

n mn  m 


(2-8) 


This  solution  may  be  exact  or  approximate,  depending  upon  the  choice  of 

e and  w . The  particular  choice  w * e is  known  as  Galerkin's  method  [9], 
n m n n 
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T-i  < 1 -1  i -A WIAr  ~ iT 


variations,  respectively.  It  is  assumed  that  the  amplitude  and 

phase  variations  of  the  excitation  are  each  normally  distributed. 

Since  magnitude  does  not  have  negative  values,  strictly  apeaking, 

the  probability  density  function  of  A,  is  a truncated  normal  dis- 

— 1 

tribution,  i.a. 


fA  U ) = — 

-i  1 J 2Ti  O 


a '(ai"an  ,2/20a 
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a^  0 


a^  < 0 


(2-14) 


where  a is  the  mean  amplitude  (a  > 0) , o is  the  standard  devi- 
mi  mi  ~ ai 

ation  of  a . K is  a constant  multiplier  which  is  obtained  by 
ai 

norm^U on  ; 


J~  V4i) 

-i 


The  error  function  erf(x)  is  defined  as: 


(2-15) 


erf (x) 


2/o 

e y dy 


(2-16) 


ai  am 

T+erf(^) 

ai 


(2-17) 


Since  phase  angle  0 is  a periodic  function  of  period  2tt,  the  probability 
density  function  of  j).  is  again  a truncated  normal  distribution. 
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= 0 


elsewhere 


where  0 is  the  mean  phase  angle,  oa  is  the  standard  deviation 
mi  0i 

of  0..  K is  a constant  multiplier.  Without  loss  of  generality, 
1 6i 

let  0 be  equal  to  zero.  Thus  the  density  function  is 
m. 


K 


0. 
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= 0 


elsewhere 


The  constant  multiplier  KQ  can  be  obtained  also  by  normilir^tiort : 
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(2-20) 
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Thus 
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(2-21) 


It  may  be  assumed  that  the  joint  density  function  of  (a^»0^)  is  a 
bivariate  normal  distribution.  , .2 


fAi0i(ai*0i)= 


K K. 
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where  p is  the  correlation  coefficient.  If  A.  and  0.  are  sta- 

~i  i 

tistically  independent,  they  are  uncorrelated,  i.e.  p = 0. 


fA.5 /VV  ‘ fA<*l>'fl,(V 
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(2-23) 


Given  a complex  random  voltage  excitation,  V^,  it  is  easy  to 
find  the  expected  value  of  if  the  amplitude  and  phase  variations 
are  assumed,  to  be  statistically  independent. 


tr 
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EP^]  = ■ da±  j d0i  fA  (ai>f  q (0±)  [a±  cos  0i  + jai  sin  0±]  (2-25) 
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E[V  ] = E[A.]E[cos  0J  + jE[A.]E[sin  6^ 


(2-26) 

(2-27) 
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(2-28) 


If  a is  small  compared  with  a , say  0 < 0.5  a for  instance, 

ai  mi  ai  mi 

the  expected  value  can  be  approximated  by 


E[A, ) = K a 
-iJ  a.  m± 


(2-29) 
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(2-31) 


If  aQ  is  small  enough,  such  that  fQ  (0.)  2 0 for  0.  = ±i j , the 
9i  1 1 


expected  value  is 


E[cos  0 ] = K~  e 
— i u . 
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ei 
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Since  sin  0 is  an  odd  function,  fQ  (0.)  is  an  even  function  over 

h 1 


f — tr ,*n ] . The  expected  value  of  sin  0^  is  zero. 


E[sin  _0  ] = 0 


(2-33) 


Thus  Eq.  (2-27)  is 


-<V2> 


E[V  ] = (K  a )(K  e 1 ) 
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(2-34) 


If  o =0  (no  amplitude  variation) , then 
ai  2 

-<°e,/2> 

E'V  ' 


(2-35) 


If  o.  =0  (no  phase  variation) , then 
6i 


E[V . ] = K a . 
— i a . m . 

i l 


(2-36) 


2.2.2  Statistics  for  Multiple  Excitations 


Let  [V  ] be  a column  matrix  with  n-elements  v.  , v„,  ...  v . If 
n 1’  2’  n 

each  element  v is  a complex  random  voltage  as  defined  in  Eq.  (2-12), 

[V  ] is  an  n-dimensional  complex  random  vector,  denoted  as  [V  ].  Or  it 
n — n 

can  be  redefined  as  a 2-n  dimensional  real  random  vector  [V»  ],  with 

— 2n 

real  random  variables,  a^ ,6^ , a_2»— 2’  as  elements* 

It  is  assumed  in  Section  2.2.1  that  the  joint-density  function 
of  is  a bivariate  normal  distribution  (Eq.  2-22).  For  simpli- 

city, it  is  assumed  that  the  random  variables  v, , v„,  ....  v are  also 

—12  — n 

statistically  independent.  It  is  possible  to  partially  isolate  each 
input  channel  of  the  array  system.  Thus,  the  random  errors  in  each 
input  excitation  can  be  assumed  statistically  independent  with  some 
justification.  Thus,  the  joint-probability  density  function  is 


fn(-l’  -2*  fV  'vl)  fV„(v2)  •••  fV  (vn) 

—1  '—Z  — n 


(2-37) 


Or  in  terms  of  a.,  0.: 
l i 

f2»'Vi2  ••••  *— 1 ‘—2  •••  V ■ fA1<al)fA2(a2>  - ^ 

If  the  random  excitations  v, ,v„,  ...  v are  not  statistically 

—1  —2  — n 

independent,  then  the  joint  probability  density  function  is  quite 

complicated.  However,  if  the  phase  variations  are  not  present,  the 

joint  density  function  v^,  v^*  •••>  is  an  n-dimensional  joint-normal 

distribution.  Let  [B  ) be  a column  matrix  with  n elements  b,  , b„,  ....  b . 

n 1 2 n 

Each  b.  is  defined  as  the  expected  value  of  v.  (see  Eq.  (2-36)).  Define 
x — 1 


i 


*t : 
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Che  covariance  matrix  [M  1 as: 

n 

[Ml  = E[  ClV  ] - [B  ])([V  ] - [B  ])T]  (2-39) 

n — n n — n n 

where  T means  transpose.  The  n-dimensional  jointly-normal  distribution 
function  can  then  be  written  as 


2.3  STATISTICS  OF  FAR  FIELD  BEAM  PATTERN 

The  calculation  of  an  electromagnetic  field  problem  can  be  simpli- 
fied to  a matrix  equation (Eq.  (2-8))  by  the  method  of  moments.  Suppose 
a matrix  [F]  of  order  mxn  is  obtained  for  a planar  array,  so  that  given 
a column  matrix  [v]  of  order  n which  represents  the  voltage  excitation 
on  each  element  of  the  array,  a column  matrix  [e]  of  order  m can  be 
obtained  which  represents  the  far  field  electric  field  intensity  at  m 
different  positions  in  space,  provided  suitable  sets  of  expansion 
and  weighting  functions  are  used.  The  matrix  equation  is 

[e]  = [F ] [ v 1 . (2-41) 

Now  consider  the  statistical  problem  of  the  planar  array. 

Eq.  (2-41)  is  still  applicable,  except  that  both  [v]  and  [e]  are  random 
vectors.  Thus, 

[e]  = [F] [v]  (2-42) 

The  elements  of  matrix  [F]  are  uniquely  determined  by  the  given  geometry 
of  the  array,  and  the  expansion  and  weighting  functions.  Thus,  each 
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far  field  quantity  e.  is  only  a linear  combination  of  random  variables 


v,  s. 

—i 


e . = l f v. 

-1  jfi  ij  rJ 


l — 1 , 2 , • * ■ | in 


(2-43) 


or,  in  another  form. 


e.  = 7 f,.a.  cos  6.  + j J f..a.  sin  0. 
“I  j,!  ijT  “3  iJ-3  “J 


(2-44) 


The  expected  value  of  the  far  field  quantities  can  easily  be 
obtained  if  the  excitations  v ^ , •••  are  assumed  statistically  inde- 


pendent. 


Et-i!  ■ v(v 


j=l 


i — 1,  2,  ■ . . , m 


(2-45) 


where  E[v.]  is  shown  in  Eq.  (2-34). 

~2 

Suppose  that  the  statistics  of  the  excitations  v ^ ...  are 

identical,  i.e.  the  probability  distribution  functions  f (v.),f  (v  ),., 

-1  1 -2  L 

are  all  identical.  Then 

n 

E[e  ] = E(v]  l f 

j = l 1 


i 1,2,  . . . , m 


(2-46) 


The  ratio  of  average  far  field  patterns  with  random  error  excitations 
to  the  error-free  patterns  (uniformly  excited)  is  obtained  by 


E[e  ] E[v]  Z f E[v] 

...  ss  — ■ — — 


e . 
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v £ f .. 

ij 


(2-47) 


If  the  v^,  v^»  •••  are  not  statistically  independent,  the  sta- 
tistics of  far  field  quantity  e.  are  complicated.  However,  if  no  phase 


variation  is  present,  then  v , v_2»  •••  are  normally  distributed  in 
r 23 1 

n-dimensions  . The  expected  value  and  the  V3.ri3.nce  of 

e^  are,  respectively, 
n 


E[eJ 


£ 

3-1 


V'V 
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and 


n n 

I l 

J-l  k=l 


f ijfikE^Vj'1 


i = 1, 
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....  ..  . .. 


1.  THEORY  OF  THE  PROGRAM 


3.1  INTRODUCTION 

Computer  program  PLAN  is  developed  for  the  purpose  of  analyzing 
the  effects  of  phase  and  amplitude  errors  on  the  far  field  beam  pattern 
of  planar  arrays  of  dipoles.  Fig.  1 shows  a typical  planar  array  of 
dipoles,  lying  in  the  x-z  plane,  with  uniform  spacing  S^,  Sz  in  the 
x,z  direction,  and  with  M2,  M^  dipoles  in  the  x,z  directions,  respec- 
tively. The  dipoles  are  centerfed  with  ideal  voltage  sources.  It  is 
assumed  that  the  statistical  variation  of  the  amplitude  and  phase  of 
the  voltage  sources  is  known.  Computer  program  PLAN  is  used  to  com- 
pute the  statistics  of  the  far  field  beam  pattern,  given  the  statistics 
of  the  voltage  sources.  In  other  words,  the  program  is  used  to  generate 
the  transfer  function  between  source  statistics  and  beam  pattern  sta- 
tistics for  the  array. 

It  is  assumed  that  the  amplitude  and  phase  of  the  signal  are 
* 

individually  Gaussian  . A nominal  phase  and  amplitude  are  specified, 
together  with  the  standard  deviation  for  amplitude  (<?A)  and  for  phase 

(ap). 

In  the  computer  program,  the  voltage  sources  of  Fig. 1 are 

selected  by  a random  number  generator  using  the  given  standard  deviations 

a , O . The  array  is  filled  in  this  manner  and  the  far  field  beam 
A r 

*It  should  be  noted  that  the  complex  voltage  is  not  Gaussian.  For  a 
complex  signal  to  be  Gaussian,  the  phase  must  be  uniformly  distributed 
(0  - 2t r). 
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1 


1 

I 

pattern  computed.  The  process  is  then  repeated  a number  of  times  1 

(typically  4Q  or  60  times*)  and  the  statistics  of  the  magnitude  and  phase 
of  the  resulting  far  field  beam  patterns  and  gains  are  computed.  This 
analysis  can  be  carried  out  at  fundamental  frequency  for  a halfwave 
dipole  (L  **  0.5A)  and  at  second  (L  = 1.0A)  and  third  (L  = 1.5A)  harmonic 
and  thus  harmonic  levels  can  be  established  relative  to  fundamental. 

Mean  beam  patterns  are  computed,  as  well  as  the  individual  beam 
patterns  and  the  expected  range  of  variation  of  the  patterns  at 
fundamental  and  harmonic  frequencies. 

Once  the  array  is  filled  up  by  the  random  number  generator, 
the  problem  is  treated  as  a deterministic  one.  The  method  of  moments 
[3,9]  is  used  to  compute  the  far  field  beam  patterns.  In  the  method 
of  moments,  the  dipole  antennas  are  subsectioned  as  shown  in  Fig.  2. 

Antenna  currents  are  not  assumed  but  are  computed.  All  mutuals  are 
taken  fully  into  account.  The  relationship  between  currents  and 
voltages  is  given  by 

[v]  = [z]]i]  O-D 

where  [z]  is  the  generalized  impedance  matrix,  z is  the  voltage 
(proportional  to  the  tangential  electric  field)  at  subsection  i, 
due  to  a current  of  unit  magnitude  flowing  on  subsection  j.  [v]  is 
a column  vector  representing  the  known  voltage  excitation  of  the  sub- 
sections and  [i]  is  a column  vector  representing  the  unknown  currents. 

* 

This  repetitive  procedure  can  be  implemented  efficiently  since  the 
second  and  subsequent  iterations  are  much  faster  than  the  first. 
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M,  = # OF  EXPANSION  FUNCTIONS  PER  DIPOLE 
M2  M3  * # OF  ANTENNAS 

M,M2M3*  # OF  UNKNOWNS  (#  OF  EXPANSION 

FUNCTIONS) 


Fig.  2.  Thin  wire  model  of  a planar  array  with  spacing  S S 
in  the  X,  Z directions.  x z 
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F.q.  (3-1)  Is  solved  by  matrix  inversion  to  obtain 

[i]  = [z]  X[v]  - [y ] [v ] (3-2) 

Thus,  the  currents  are  found  by  matrix  inversion;  the  far  fields  are 
then  determined  from  the  currents. 

There  are  many  current  representations  possible  within  the  method 
of  moments.  In  program  PLAN,  a Galerkin  [9]  solution  with  piecewise 
sinusoidal  expansion  and  weighting  functions  is  used.  This  yields  a 
very  accurate  solution  and  permits  one  to  use  a relatively  small  number 
of  current  unknowns. 

For  the  planar  array  of  fubsectioned  dipoles  shown  in  Fig.  2,  the 
generalized  impedance  matrix  has  certain  special  features.  There  are 
a number  of  symmetries  present  and  the  impedance  matrix  has  the  form  of 
a block-Toeplitz  matrix.  In  addition,  the  voltage  excitation  column 
matrix  [v]  consists  mostly  of  zeros  since  most  subsections  are  not  ex- 
cited. Subroutines  are  developed  to  take  advantage  of  these  features. 

The  result  is  an  efficient  computer  program  which  may  be  used  to  treat 
up  to  100  dipoles  at  the  third  harmonic  (L  = 1.5AA.  The  analysis  can 
be  carried  out  at  any  frequency. 

The  program  PLAN  can  also  be  used  for  deterministic  array  analy- 
sis, one  need  merely  set  = Op  = 0 and  just  carry  out  one  computation 
(no  iterations).  The  program  is  reasonably  efficient  when  used  in  this 
manner. 

Statistical  effects  in  antennas  have  been  treated  previously  [4-7] 
but  the  previous  analyses  have  not  taken  the  effect  of  mutuals  fully 
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into  account.  The  program  PLAN  permits  one  to  include  mutuals  or  not. 
Results  show,  as  expected,  that  the  effect  of  mutuals  is  significant. 

If  mutuals  are  excluded,  the  array  problem  becomes  quite  a simple  one. 

In  this  introductory  section,  the  basic  methods  employed  in  pro- 
gram PLAN  are  outlined  briefly.  The  details  of  these  methods  are  covered 
in  sections  3.2  to  3.7. 

3.2  GENERATING  RANDOM  ERRORS  IN  THE  EXCITATIONS 

The  process  of  filling  up  the  array  of  Fig.  1 with  voltage  sources 
chosen  on  a statistical  basis  was  briefly  described  in  section  3.1.  It 
is  assumed  that  the  magnitude  and  phase  of  the  voltage  excitation  v are 
both  normally  distributed  (truncated),  e.g. 

V = Ae1-  (3-3) 

See  Section  2,2  for  details  of  the  statistical  formulation. 

A subroutine  is  used  to  generate  normally  distributed  random 
numbers  for  amplitude  and  phase  by  the  power  residue  method  [10].  Nega- 
tive random  numbers  generated  for  A as  well  as  the  random  numbers  out- 
side the  [0  -7i,  0 +tt]  range  for  I?; are  discarded, 
mm 

A and  _0  are  also  assumed  to  be  statistically  independent  in  this 
program.  Thus,  two  normally  distributed  random  numbers  a and  0 are 
generated  separately  by  the  subroutine,  and  combined  to  form  the  voltage 
excitation  v (Eq.  3-3). 

In  fact,  a bivariate  sampling  process  can  be  chosen  to  generate 
two  jointly  normal  random  numbers  with  specified  correlation  coeffi- 
cient p,  (-1  < p < 1),  which  enables  one  to  relate  the  amplitude  error 

~ ~ 
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I 
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to  the  phase  error  in  any  degree,  i.e.  from  no  correlation  to  linearly 
dependent  . 

The  bivariate  sampling  procedure  consists  simply  of  taking  two 

random  numbers  y^  and  y^,  say,  from  the  normal  distribution  with  mean 

zero  and  standard  deviation  1,  and  forming 

a = a + ay  (3-4a) 

m a 1 

0=0  + pa  + s 1-p2  a.y  (3-4b) 

m u z 

The  random  number  generator  routine  in  this  program  is  a modifi- 

★ 

cation  of  one  appearing  in  the  IBM  Scientific  Subroutine  Package  [11]. 

The  array  is  thus  filled  up  with  sources  chosen  by  the  sub- 
routine. The  analysis  is  carried  out  to  obtain  the  far-field  beam 
pattern.  Then  the  array  is  filled  up  again  and  the  analysis  repeated. 

A number  of  such  analysis  are  carried  out  and  the  statistics  of  the 
resulting  beam  pattern  are  then  computed. 

In  the  computer  program,  the  phase  and  amplitude  error  of  an 
individual  source  are  assumed  to  be  independent  (P  = 0).  In  addition, 
the  voltages  v ^ , v,,,  . . . , v^  are  assumed  to  be  independent.  As 
noted,  a bivariate  normal  distribution  ( p 4 0)  can  readily  be  treated, 

and  in  addition,  correlation  between  sources  v,  ...  v can  be  taken 

—1  — n 

into  account  using  the  covariance  matrix  (some  subroutines  are  already 
available  for  this  purpose  in  IMSL).  Also,  the  assumption  of  a Gaussian 
distribution  is  not  critical  to  the  analysis.  Any  other  distribution 

* 

The  newest  version  of  the  Scientific  Subroutine  Package  is  the  Inter- 
national Mathematical  Statistical  Library  (IMSL)  which  includes  more 
random  number  generating  subroutines. 
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can  be  simulated;  it  is  necessary  only  to  substitute  a new  random 
number  generator  with  the  desired  statistics. 


3.3  SETTING  UP  THE  MATRIX  EQUATION 


The  first  step  in  the  method  of  moments  technique  is  to  set  up 
the  matrix  equation  (Eq.  (3-1)),  i.e.  to  fill  the  matrix  [z].  Within 
the  method  of  moments,  there  are  many  possible  methods  for  the  compu- 
tation of  the  typical  element  z ^ of  the  generalized  impedance  matrix. 
If  the  same  functions  are  used  for  expansion  and  weighting  [9],  then 
a Galerkin  (variational)  solution  results.  In  program  PLAN,  a 
Galerkin  solution  with  piecewise  sinusoidal  expansion  and  weighting 
functions  is  used,  similar  to  those  developed  by  Richmond  [12]  and 
Sarkar  [13].  The  subroutine  used  in  PLAN  is  a modification  of  the 
program  of  Sarkar  [13],  the  principal  modification  being  the  addition 
of  a simple  method  of  filling  the  matrix  from  a knowledge  of  the  first 
row  only.  This  simplification  results  from  the  symmetries  present  in 
the  array;  the  details  of  these  symmetries  are  discussed  in  section 


The  use  of  piecewise  sinusoidal  expansion  and  weighting  functions 
yields  a very  accurate  solution  of  the  array  problem  and  permits  one 
to  use  a relatively  small  number  of  unknowns.  Three,  five  and  seven 
expansion  functions  (corresponding  to  3,5,7  unknowns)  per  dipole  are 
used  for  L - G.5,  1.0,  1.5A,  respectively. 
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A Galerkin  solution  results  in  an  impedance  matrix  which,  is 
symmetric  about  the  main  diagonal  i.e.  z^  = Zji"  leads  to  a 

savings  in  inversion  time  and  storage. 

The  piecewise  sinusoidal  function  used  can  be  written  as 


A sin  k(Z  - Z ) 

rn(Z)  “ WZ  sin  k(Z  7) 

n n-1 

sin  k(Z  _ - Z) 

~ n+1 

UZ  sin  k(Z  - Z ) 
n+1  n 


Z < Z < Z (3-5a) 

n-1  — n 


Z < Z < Z 
n — n+1 


<3-5b) 


where  p is  a unit  vector  in  the  z direction,  and  k = 2ir/X.  Each 
expansion  function  exists  over  two  adjacent  subsections;  the  re- 
sulting current  approximation  for  a dipole  is  illustrated  in  Fig.  3 • 
Note  that  the  number  of  subsections  in  this  case  is  always  one  more 
than  the  number  of  expansion  functions. 


3.4  THE  GENERALIZED  IMPEDANCE  MATRIX  OF  A PLANAR  ARRAY 

3.4.1  Impedance  Matrix  of  an  Isolated  Dipole 

Fig.  4a  shows  an  isolated  dipole  antenna  and  Fig.  4b  shows  the 

corresponding  filamentary  model  commonly  used  in  the  method  of  moments 

* 

for  thin-wire  antennas  . The  number  of  subsections  is  n+1;  the  number 
of  expansion  functions  is  n.  The  generalized  impedance  matrix  for 
this  antenna  may  be  written  as 

r*, 

[z]  = 


11 


In 


&nl 


nn  J 


(3-6) 


Note  that  the  dipole  radius  a is  taken  fully  into  account  by  applying 
boundary  conditions  at  the  wire  surface. 
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AL  = SUBSECTION  LENGTH 


♦ OF  SUBSECTIONS  * M,  + I 

# OF  EXPANSION  FUNCTIONS  * M,=  n 

2 j j = z j i (FOR  A GALERKIN  SOLUTION) 

2 i j = 2(i  + k)(j  + k ) (TOEPLITZ  PROPERTY) 


Fig.  4.  An  isolated  dipole  (a)  dipole  (b)  filamentary 
model  (c)  expansion  functions. 
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For  the  dipole  antenna  shown  in  Fig.  4,  z^  « z^,  *u  = z etc.,  or, 
in  general 

Zrs  Z(r+k)(s+k)  C3-7» 

where  k is  an  integer  and  r+k  n,  s+k  n.  This  property  (Eq.  C3-7)) 
defines  the  matrix  as  a Toeplitz  matrix  [14-18).  In  addition,  since  a 
Galerkin  solution  is  used,  the  matrix  [z]  is  symmetric  about  the  main 
diagonal : 


z - z 
rs  sr 


* 1 £ r ,s  <_  n 


(3-8) 


The  result  is  a symmetric  Toeplitz  matrix;  the  elements  of  the  matrix 
are  repeated  along  diagonal  lines: 


[*]  - 


Z0  Z1 


Z1  zo 


' z ,z 
n - 1 n- 2 


n-i 


(3-9) 


where  z defines  the  self-impedance  term  of  one  subsection  and  z, 

0 I r-s | 

defines  the  mutual  impedance  term  between  two  subsections  r and  s. 

Note  that  the  symmetric  Toeplitz  matrix  shown  above  can  be 
generated  by  the  first  row  or  column.  This  greatly  simplifies  the 
fill  time  (fill  time  is  proportional  to  n rather  than  n ) and  reduces 
the  computer  storage  required.  A non-symmetric  Toeplitz  matrix 
would  be  generated  by  knowledge  of  the  first  row  and  the  first 
column. 

Toeplitz  matrices  occur  in  the  theory  of  discrete  random  processes. 


r; 


Several  special  algorithms  have  been  developed  for  the  inversion  of 
such  matrices  [15-16],  The  matrix  [y]  is  defined  as  the  inverse  of 


[y]  - [z]  1 = I 


The  Toeplitz  symmetry  of  the  [z]  matrix  is  lost  upon  inversion, 
but  there  is  a fourfold  symmetry  remaining.  This  can  be  explained  by 
the  definition  of  y and  Fig.  5. 

y^  = — (with  all  other  v’s  equal  to  zero).  This  represents 
the  input  admittance  at  port  i (defined  in  terms  of  the  ith  expansion 
function)  with  all  other  ports  short  circuited. 


y = — (with  all  other  v's  equal  to  zero).  This  represents 

J VJ 

the  transfer  admittance  between  ports  i and  j with  all  other  ports 
short  circuited. 

[y]  is  symmetrical  about  the  main  diagonal  because  a Galerkin 
solution  is  used,  i.e.  z^  ■ z imPl*es 

y±s " yji  <3'n> 

In  addition,  [y]  is  symmetrical  about  the  secondary  diagonal 
(the  diagonal  from  y^n  to  y^) : 

^ij  ^(n-j+1)  (n-i+1)  for  (l<^i<^n,  1 £ j <_  n).(3-12) 


This  property  can  be  expTained  by  reference  to  Fig.  5;  it  is  a con- 
sequence of  physical  symmetry  about  the  center  of  the  dipole. 


Thus  the  [y]  has  a fourfold  symmetry  (about  main  and  secondary  diagonals) 
as  illustrated  in  Fig.  5* 

In  addition  to  the  symmetries  in  [z]  and  [y] , which  may  be 
called  redundancies  in  impedance  and  admittance  elements,  and  which 
are  independent  of  excitation,  there  may  be  symmetries  in  the  currents. 
For  instance,  if  the  dipole  is  center  fed,  there  is  an  additional 
two-fold  symmetry  in  currents.  Current  symmetries  can  readily  be 
taken  into  account  by  using  special  symmetry  subroutines  which  reduce 
the  number  of  unknowns  if  the  current  symmetries  present  are  specified. 

In  program  PLAN  there  are  no  such  symmetries  present  because  of  the 
mutual  effects  and  because  of  the  randomness  of  the  voltage  excita- 
tions. 

3.4.2  Impedance  Matrix  of  a Linear  Array 

Fig.  6a  shows  an  equally  spaced  linear  array  with  N+l  identical 
dipoles.  The  corresponding  thin-wire  model  is  shown  in  Fig.  6b  . Ml 
expansion  functions  are  used  for  each  dipole. 

The  generalized  impedance  matrix  is  a square  matrix  of  order 
Mix  (N+l).  For  a suitable  ordering  of  expansion  functions  , the 
generalized  impedance  matrix  can  be  partioned  as 

*11  *12  *1( N+l) 

*21  *22 

Iz]  " : : 0-13) 

• • 

*N+1)  1 *N+lxN+l> 

L J 

* 

The  expansion  functions  are  ordered  continuously  over  each  dipole, 
beginning  at  a corresponding  point  on  each  dipole  as  shown  in  Fig. 6b. 
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where  each  submatrix  z is  a square  matrix  of  order  Ml.  represents 

interaction  among  portions  of  antenna  i and  z represents  interactions 

A 

between  antennas  i and  j . z^  is  thus  a matrix  identical  to  that  of 
Eq.  (3-6)  for  the  isolated  dipole.  Since  all  antennas  are  identical, 
all  are  identical,  i.e. 


A 

*ii  = *0 


for  i = 1,2,. .. ,N+1 


(3-14 


where  = designates  a definition.  All  z are  thus  symmetric  Toeplitz 
matrices  (symmetric  because  a Galerkin  solution  is  used) . 

The  typical  submatrix  z^  involves  the  interactions  between 
antennas  i and  j.  The  matrix  [z]  is  symmetric  about  the  main  diagonal 
and  thus 


*ij  = % 


(3-15) 


in  general.  For  the  particular  case  in  which  the  line  of  the  array 
is  perpendicular  to  the  dipoles  (the  "perpendicular  line"  case  of 


Fig.  6a) 


ill  ' ’ Ail 


(3-16) 


where  superscript  T indicates  the  transpose  operation. 


For  members  i,j  of  the  linear  arrays  of  Fig.6a,  z depends 
only  on  the  distance  between  antenna  centers  and  thus 


Note  that  each  generalized  impedance  element  represents  the  inter- 
action between  two  current  expansions.  For  z..,  both  expansion  func- 
tions are  located  on  the  same  antenna  and,  for  z. . , the  expansion 
functions  are  located  on  different  antennas  (i^j)^ 
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■■■■■■■■SMI 


‘ij  “ *|i-j| 


for  (1  < i < N+l,  1 < j < N+l) , 


0-17) 


z^.  Is  itself  a Toeplitz  matrix,  since  the  elements  of  the  matrix 


depend  only  on  the  difference  between  subscripts  t 
is  itself  a symmetric  Toeplitz  matrix  for  the  "perpendicular  line" 
case. 

The  resulting  entire  impedance  matrix  [z]  is  then 
z. 


*0 


• • 


% 


T 

I *l 

-0  * * * %-l 

[z]  = 

• 

• 

• 

1 T 

• 

T 

i -?N 

-N-l  ’ ’ ~0  . 

simplified  in 

the  "perpendicular 

*0 

h • ’ ~N 

1 

1 ~1 

^0  ~N-1 

Z = 

• 

~N-1  ~0 

for  Eq.  (3-19 

' , all  submatrices 

(3-18) 


(3-19) 


matrices. 

A matrix  of  the  form  (3-19)  is  called  a "block-Toeplitz"  matrix 
since  Eq.  (3-19)  is  of  the  same  form  as  Eq.  (3-9)  if  elements  are 
replaced  with  submatrices,  i.e.  the  blocks  are  repeated  along  diagonal 


lines. 


The  entire  matrix  of  (3-18)  can  be  generated  by  the  first  row 
and  the  first  column.  The  entire  matrix  of  (3-19)  can  be  generated 
by  the  first  row  or  the  first  column.  This  latter  statement  can  be 
verified  by  noting  that  the  matrix  can  be  generated  by  the  first 
row  of  (N+l)  submatrices  (z^,  z^,  ...  z^) . Each  submatrix  can  be 
generated  from  its  first  row  since  it  is  individually  Toeplitz. 

Thus  the  number  of  independent  impedance  elements  is  equal  to  the 
number  of  elements  in  the  first  row,  which  is  equal  co  Ml (tn-1)  or 


or  , 


[y] , the  inverse  of  [z],  is  designated  as  follows: 


r?oo 

y01 

y02  ' ' ’ 

y0N 

yio 

yll 

y12 

ZlN 

j 

[ yN0 

yNl 

yN2 

yNN 

(3-20) 


Since  a Galerkin  solution  is  used,  the  full  matrix  is  symmetrical  about 
the  main  diagonal,  and  the  blocks  y are  symmetrical  about  the  main 
diagonal  after  transposition,  i.e. 

T 

«•  / /-v  . . . ^ i \ 


" y^i 


for  (0  <_  i,  j £ N) 


In  addition  there  is  a fourfold  symmetry  of  the  blocks: 


^ij  = ? (N-i)  (N- j ) “ y (N-j ) (N-i) 


(3-21) 


(3-22  ) 


Thus  the  entire  y matrix  can  be  generated  from  one-quarter  of  the 
blocks  (see  Fig.  7).  The  above  is  valid  for  both  the  "perpendicular 
line"  and  "skewed  line"  cases. 
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Fig.  7.  An  equally  spaced  linear  array  with  N+l  identical  dipoles  - 
symmetries  of  the  [Y]  matrix  (M„  = 5). 


The  individual  blocks  also  have  some  redundancies  of  admittance 


elements.  Consider  first  the  blocks  y . These  blocks  have  a fourfold 
symmetry  about  the  main  diagonal  for  the  "perpendicular  line"  (see 
Fig.  6a).  por  the  "skewed  line"  case,  there  is  a twofold  symmetry 
about  the  main  diagonal  (see  Fig.  6c). 

Next  consider  the  blocks  y . First  of  all  the  blocks  along  the 
secondary  diagonal  have  a fourfold  symmetry  in  the  'perpendicular  line" 
case  which  is  lost  in  the  "skewed  line"  case.  The  only  blocks  re- 
maining are  those  both  off  the  main  and  secondary  diagonals.  These 
blocks  have  a different  type  of  symmetry;  namely  a twofold  reflection 
symmetry  about  the  center  element  in  the  "perpendicular  line"  case 
(see  Fig.  7).  This  symmetry  is  lost  in  the  "skewed  line"  case. 

In  summary,  there  are  many  redundancies  in  the  [z]  and  [y] 
matrices.  For  the  "perpendicular  line"  case  [z]  can  be  generated 
by  the  first  row;  in  other  words  there  are  M^(N+1)  independent 
elements  of  [ z ] . For  the  "perpendicular  line"  case,  there  is  a 
fourfold  symmetry  of  blocks,  and  at  least  a twofold  symmetry  in  the 
blocks  themselves.  [y]  can  thus  be  generated  by  somewhat  less  than 
one  eighth  of  the  elements. 


3 .4.3  Impedance  Matrix  of  a Rectangular  Planar  Array. 

The  impedance  and  admittance  matrices  of  a planar  array  are, 
of  course,  considerably  more  complex  than  those  of  the  linear  array. 
It  would  require  considerable  space  to  present  all  of  the  relation- 
ships in  detail.  Some  of  the  essential  relationships  can  be  outlined 


briefly  by  considering  the  planar  array  as  a linear  array  of  elements, 
each  element  being  a linear  array  itself. 


For  example  the  rectangular  planar  array  of  Fig.  8 can  be  treated 
as  a linear  array  by  taking  either  (a)  a row  of  dipoles  or  (b)  a 
column  of  dipoles  as  an  element.  The  resulting  linear  arrays  are 
shown  in  Figs,  9a,  9b. 

The  impedance  matrix  of  the  array  may  be  written  as  follows: 


-=1(N+1) 

Z2(N+-1) 


(3-23> 


^CN+1)1  Z(N+1)2 


:(N+1)(N+1) 


where  Z denotes  interaction  among  expansion  functions  of  array  i 
and  denotes  interactions  between  arrays  i and  j . 

Since  all  arrays  are  identical  Z^  may  be  defined  as  Zq  and 
Z^  may  be  defined  as  Z|^  ^ | , in  analogy  to  Eqs.  (3-14)  and  (3-17). 
The  result  is  a matrix  of  the  block- Toeplitz  form: 


zo 

?1 

?2  •* 

% 

^1 

^0 

^1 

~N-1 

024) 


\h  Vl 


~N-2 •'  V 


Note  that  Zq  is  the  impedance  matrix  of  an  array  and  it  is  thus  also 
of  the  block-Toeplitz  form.  The  same  is  true  of  the  off-diagonal 
blocks.  Since  each  block  of  (3-24)  is  block-Toeplitz,  the  impedance 
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matrix  of  a planar  array  may  be  termed  a block-block  Toeplitz  matrix. 


If  the  planar  array  is  square  (M2  ■ M3)  then  the  matrix  inversion 
problem  can  be  treated  equally  well  by  taking  (a)  rows  or  (b)  columns 
of  dipoles  as  elements.  If  the  planar  array  is  not  a square,  the  shorter 
side  of  the  rectangular  array  should  be  taken  as  an  element.  This  gives 
the  best  results  as  explained  in  Appendix  A. 

For  example,  a planar  array  of  3 by  4 dipoles  is  sh^wn  in  Fig.  9. 

It  is  preferable  to  choose  the  3-dipole  column  as  an  element  rather  than 
the  4-dipole  row.  The  resulting  impedance  matrix  will  contain  16  sub- 
matrices in  comparison  with  the  9 submatrices  resulting  from  the  choice 
of  the  4-dipole  row  as  an  element. 

Note  that  by  analogy  to  the  previous  section  the  generalized  impe- 
dance matrix  of  the  rectangular  planar  array  can  be  generated  from  the 
first  row  or  column  of  [ z ] . 

The  admittance  matrix  [y]  has  a fourfold  symmetry  about  the  main 
and  secondary  diagonals,  as  before,  for  the  rectangular  planar  array. 


3.5  SOLUTION  OF  THE  MATRIX  EQUATION 

The  matrix  equation  for  the  planar  array, 

[v]  = [z][i]  (3-2 
is  first  set  up  using  the  methods  described  in  section  3.3  and  the  sym- 
metries outlined  in  3.4.  The  solution 

[i]  = [z]  1[v]  = [y] [v]  (3-2 
is  obtained  in  two  steps.  The  first  step  involves  the  inversion  of  the 
generalized  impedance  matrix  [z]  which  is  of  the  block-block-Toeplitz 
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torm.  This  is  accomplished  by  a special  routine  which  is  described 
in  detail  in  Appendix  A.  Essentially  the  routine  is  derived  from  the 
method  used  by  Sinnott  for  linear  arrays.  In  the  routine  described 
in  Appendix  A,  there  are  three  additions  to  the  method  of  Sinnott: 

(1)  The  method  is  applied  to  planar  arrays  by  considering 
either  a column  or  row  of  elements  as  an  element.  Details  are  given 
in  Appendix  A. 

(2)  The  zeros  of  the  excitation  matrix  are  utilized  to  simplify 
the  inversion  and  reduce  the  core  storage.  Most  of  the  elements  of 
the  voltage  excitation  matrix  [v]  are  zeros.  In  the  final  solution 
(3-26)  then,  only  the  columns  of  the  matrix  [y]  corresponding  to  the 
non-zero  elements  of  [v]  are  required.  In  other  words,  the  entire 

[y]  matrix  is  not  required.  The  inversion  routine  can  be  correspond- 
ingly simplified.  In  addition,  the  matrix  multiplication  (3-26) 
can  be  simplified  by  omitting  the  corresponding  multiplications.  This 
is  explained  further  in  Appendix  A. 

(3)  The  main  (core)  computer  memory  requirements  are  reduced 
by  the  use  of  peripheral  devices  (disk  storage).  Two  disc  files  are 
used  as  temporary  storage.  Details  are  given  in  Appendix  A. 
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3.6  COMPUTATION  OF  CURRENTS  - ZEROS  OF  THE  EXCITATION  MATRIX 


As  noted  in  section  3.5,  most  of  the  elements  of  the  voltage 
excitation  matrix  [v]  are  zeros.  Therefore,  only  a few  columns  of 
[y]  are  required  in  the  matrix  multiplication  [y][v],  namely  those 
corresponding  to  the  non-zero  elements  of  [v]. 

To  illustrate  the  method,  consider  [y]  and  [v]  matrices  where 
[y]  is  an  nxn  matrix  and  [v]  is  an  n-element  column  vector.  There 
are  only  certain  non-zero  elements  in  [v].  Suppose  that  there  are  k 
such  non-zero  elements.  The  matrix  equation  0-26)  may  then  be  written 


ID  - [YjTj  ...  YkJ  Kl 


0-27  ) 


where  the  elements  v^  ...  v^  are  non-zero.  For  simplicity,  use  a 
numbering  system  for  expansion  function^ .unknown  currents  and  voltages 
which  starts  with  the  driven  sections  1 ...  k.  Then  the  matrix  equation 
becomes 


fi]  - [Y^  ...  Yr  ...  YnJ  r V]L1 


(•3-28> 


where  Y^  is  the  ith  column  of  the  [y]  matrix.  Eq.  (3-26)  may,  of  course, 
be  simplified  to 


i 


4) 


r — -i  — I*  ■ - 


.•mgav  SH 


[i]  = [Y1  Y2  ...  Yk]  rVl* 


(3-29) 


Only  the  columns  Y^  ...  Y^  of  the  [y],  matrix  need  be  computed  and  only 
the  multiplications  indicated  in  Eq.  (3-29)  need  be  carried  out. 

Thus  the  inversion  routine  and  the  computation  of  currents  can  be 
simplified  by  taking  advantage  of  the  zeros  of  the  excitation  matrix. 

This  is  especially  useful  for  cases  where  the  number  of  expansion  func- 
tions per  dipole  is  large,  for  instance,  for  the  third  harmonic  analysis 
of  a dipole  which  is  one-half  wavelength  long  at  the  fundamental  frequency. 
In  the  computed  results,  seven  subsections  per  dipole  have  typically  been 
used  for  the  third  harnonic  analysis.  For  this  choice,  only  one  seventh 
of  the  columns  of  the  [y]  matrix  need  be  computed  or  stored  and  only 
one  seventh  of  the  multiplications  of  Eq.  (3-26)  need  be  carried  out. 


3.7  COMPUTATION  OF  FAR  FIELDS 

Once  the  currents  [i]  have  been  obtained,  the  far  fields  may  be 
found  directly  by  an  additional  matrix  multiplication.  For  z-directed 


dipoles 


[E]  = [F] [i] 


f Jt 
» '• 

(3-30) 


where  [E]  is  an  n-element  column  vector  representing  values  of  the 


electric  field  at  n different  points  in  the  far  field.  F_  is  the  elec- 
tric fieJd  at  point  i (in  the  far  field)  due  to  current  expansion  func- 
tion j of  unit  magnitude.  Thus,  f F]  is  proportioned  to  the  far  field 
form  of  the  generalized  impedance  matrix  [z]. 


- — .•  _ | 


42 


Once  the  far  field  [E]  has  been  calculated,  the  array  gain,  G, 
may  also  be  found  by 


A Air 


.2,  ,2 


where  P^n  is  the  total  input  power  to  the  array,  which  can  be  obtained 
simply  by  adding  the  input  powers  to  each  element  of  the  array. 

N+l 

P.  = ReY  V.I.  (3- 

m l l 

i=l 


There  are  many  subroutines  available  for  the  computation  of  far 
fields.  The  subroutine  developed  by  Sarkar  [13]  has  been  utilized 
here.  This  subroutine  computes  the  far  fields  due  to  a piecewise  sinu- 
soidal current  distribution. 


3.8  PROCESSING  OF  THE  FAR  FIELD  DATA 

The  sampling  procedure  for  filling  up  the  array  is  described  in 
section  3.2.  Given  the  statistics  of  the  sources,  i.e.  , the  standard 
deviations  a^,  Op  in  amplitude  and  phase,  a random  number  generator  is 
used  to  select  a set  of  sources  for  the  array.  Then  the  matrix  equation 
is  solved  and  the  far  fields  are  computed. 

Next  a new  set  of  sources  is  selected  for  the  array  and  the  whole 
process  is  repeated.  Of  course,  on  the  second  round,  [y]  need  not  be 
recomputed,  but  only  the  matrix  multiplications  [y][v]  and  [F][i]  need 
be  carried  out.  This  is  fortunate,  since  the  matrix  multiplications 
are  very  rapid  as  compared  to  the  original  matrix  inversion  [y]  = [z]  \ 
The  process  is  repeated  further.  Finally,  K sets  of  array  sources 
have  been  selected  and  K sets  of  currents  and  far  field  data  have  been 
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obtained.  Only  one  matrix  inversion  [y]  = [z]  has  been  performed; 

K.  matrix  multiplications  [y][i]  and  [F][i]  have  been  carried  out.  Typi- 
cal values  for  K are  40,  60  in  the  data  shown  in  section  4.  K is  speci- 
fied by  the  user  as  part  of  the  input  data. 

Now  the  far  field  data  must  be  processed.  First  the  complex  far 
field  data  is  converted  to  magnitude  form.  At  each  point  in  the  far  field, 
the  average  magnitude  of  the  far  field  and  gain  is  obtained,  yielding  the 
average  beam  pattern  and  gain  magnitude.  Next  the  magnitude  data  at  each 
far  field  point  is  examined  and  the  standard  deviation  a is  obtained. 

Mean  (M)  and  Mean  minus  sigma  (M-a)  beam  patterns  are  then  plotted  for 
the  array  and  the  given  source  statistics.  The  beam  pattern  M-o  reveals 
the  expected  range  of  variation  for  the  given  problem.  One  can,  of  course, 
examine  further  the  data  at  each  far  field  point  and  readily  obtain  any 
statistics  desired. 

For  example,  the  covariance  matrix  of  the  far  field  points  tells 
the  correlation  between  each  far  field  point.  The  Tchebycheff  inequal- 
ity equation,  i.e. 

2 

P(|x  * *J  > O 1 ’ (e  > 0)  (3-33) 

z 


gives  a range  of  variation  with  specified  probability  of  occurrence  within 
that  range. 

All  beam  patterns  are  normalized,  but  in  addition,  the  absolute 
level  of  the  maximum  is  indicated  by  a normalization  number  which  is 
printed  on  each  beam  pattern  plot. 

The  procedure  outlined  above  can  be  carried  out  at  harmonic  (2Fq,3Fq) 
as  well  as  fundamental  (F^)  frequencies  and  thus  the  levels  of  harmonics 
relative  to  fundamental  may  be  obtained. 


i 
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3.9  CAIN  DEGRADATION 

Random  phase  and  amplitude  errors  in  an  array  have  the  effect  of 

reducing  directive  gain.  For  large  phase  and  amplitude  errors, this  gain 

reduction  is  quite  significant.  This  effect  is  important  in  the  analysis 

of  2nd  and  3rd  harmonic  emissions  in  arrays  since  phase  and  amplitude 

errors  are  in  general  significant  at  these  frequencies  [1],  Thus,  there 

will  be  in  general  at  2nd  and  3rd  harmonics  a significant  reduction  in 

gain  over  that  obtained  from  an  optimally-excited  array  (such  as  a uniformly 

excited  array).  The  gain  of  a uniformly  excited  (non-random)  array  is 

denoted  G . 

o 

There  are  several  definitions  offcgain  which  are  of  interest  here. 

Suppose,  for  instance,  that  an  array  is  excited  for  broadside  radiation. 

Consider  a set  of  typical  sample  beam  patterns  such  as  those  shown  in 

Fig.  18-20.  One  important  quantity  would  be  the  expected  directive  gain 

for  a given  location  (9  ,0  ). 

o o 

E[G(0  ,0  )]  = Expected  directive  gain  at  9 0 (3-34) 

— n n n n 


The  expected  gain  will,  of  course,  vary  with  location.  It  represents 

.he  "expected"  interference  for  an  observer  at  (0  ,0  ). 

o o 

The  maximum  expected  gain  is  also  of  interest,  where  the  maximum 
expected  gain  is  defined  as  the  maximum  value  of  the  expected  gain  over 
all  angles  ®O>0Q*  This  quantity  is  denoted  G^: 

G^  = Maximum  Expected  Gain 

= MAX [ E [ G ( 0 ,0  )]]  (3-35) 

— oo 

(where  the  above  quantity  is  the  maximum  over  all  angles  9q,0o). 
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For  the  samples  of  Figs.  18-20,  these  values  of  maximum  directive  gain 
occur  at  different  angles.  Maximum  gain  is  thus  defined  as  the  maximum 
directive  gain  of  an  individual  sample  pattern. 

To  compute  maximum  expected  gain,  one  computes  an  expected  gain 
(by  averaging  over  the  sample  patterns)  for  each  angle  ®O»0O  and  then 
chooses  the  largest  of  these  expected  gains.  For  nominal  broadside 
excitation,  for  instance,  maximum  expected  gain  will  occur  at  broadside. 

thus  represents  the  expected  gain  at  a particular  position  at  which 
the  expected  gain  is  largest. 

Another  important  quantity  is  the  "expected  maximum  gain"  which 
is  defined  as  follows.  Consider  the  sample  patterns  such  as  those  of 
Figs.  18-20.  For  a particular  sample  pattern,  there  is  a value  of 
maximum  directive  gain.  The  expected  maximum  gain  is  then  an  average  over 
all  samples: 


= Expected  maximum  gain 

= E[G  ] 

—max 


(3-36) 


where  G is  the  maximum  directive  gain  of  a particular  sample  pattern 

and  G is  the  associated  random  variable.  The  expected  maximum  gain 
—max  s 

thus  represents  an  important  measure  of  worst-case  interference  capa- 
bility. 

The  quantities  Gq,  E [G^C&q , ) ] , G^ , G^  differ  significantly  in 
many  cases.  To  adequately  analyze  an  array  with  random  errors.it  may 
be  necessary  to  take  all  these  quantities  into  account. 


3. 10  ADDITIONAL  BEAM  PATTERN  QUANTITIES 


The  quantities  discussed  In  sections  3. 1-3.9  are  computed  in  program 


PLAN.  There  are,  of  course,  other  beam  pattern  quantities  which  could 


readily  be  computed  from  the  sample  patterns  by  modification  of  program 


PLAN.  Some  of  these  are  (1)  expected  beam  pointing  error,  (2)  expected 


beam  broadening,  (3)  expected  ratio  of  power  in  main  beam  to  power  in 


3.11  SOME  SPECIAL  PLANAR  ARRAYS 


Computer  program  PLAN  (described  in  section  3)  is  designed  to 
solve  problems  for  the  class  of  rectangular  or  square  planar  arrays. 

It  can  also  readily  be  modified  to  treat  some  other  special  planar 
arrays.  It  is  the  purpose  of  this  section  to  describe  these  special 
arrays  and  the  modifications  necessary  for  their  analysis. 

First,  consider  the  class  of  triangular  planar  arrays  (Fig.  10) 
which  is  often  encountered  in  practical  design.  If  two  adjacent 
columns  (Fig.  11a)  of  antennas  are  taken  as  an  element,  the  result 
is  a linear  array  (perpendicular  case).  If  two  adjacent  rows  (Fig.  lib) 
are  taken  as  an  element,  the  result  is  also  a linear  array  (skewed 
line  case).  In  either  case,  the  generalized  impedance  matrix  can  be 


partioned  as: 


-•1  2 ~1  (N+l) 


~(N+1)  1 ~(N+1)  (N+l) 


where  Z^  represents  the  generalized  impedance  of  element  i which 
consists  of  a double  row  or  column  linear  array.  Z^  represents  the 
interactions  between  elements  i and  j , each  of  which  is  a double  row 
(or  column)  linear  array. 

Note  that  submatrices  Z^  depend  only  on  the  magnitude  of  the 
difference  between  i and  j : 
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^ij  * ~ji  = C 

for  (L  _<  i,  j N+l)  and  [z]  may  be  expressed  in  the  following  form 
which  coincides  with  (3-24) 


h 


[z]  -j. 


-N-li 


!^N  ~N-1 


The  impedance  matrix  [z]  is -a  block-Toeplitz  matrix.  The  sub- 
matrices Z^,  0 <_  i N are  not  necessarily  block-Toeplitz  matrices, 
because  a double  row  (or  column)  element  has  been  used.  The  computer 
program  PLAN  is  already  designed  to  treat  matrices  of  the  form  of  (3-38) 

and  thus  the  only  change  required  is  to  generate  the  elements  Zn  to  Z„ 

^-U  — N 

for  the  particular  geometries  of  Figs.  11  a,  b.  There  are,  of  course, 
several  different  ways  of  generating  the  elements  Zq  to  Z^.  One  could, 
for  instance,  use  one  of  the  general  subroutines  [13,20,21]  available 
for  arbitrary  arrangements  of  elements. 

Figs.  12  and  13  show  a parallelogram  and  an  unequally  spaced  planar 
array.  These  arrays  can  also  be  treated  by  considering  single  or  double 
rows  or  columns  as  elements.  The  generalized  impedance  matrix  is  of 
the  form  (3-36)  and  one  need  only  add  a subroutine  for  the  generation 


of  the  elements  Zn  to  Z„  of  the  matrix. 
~0  —N 
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4.  PROGRAM  DESCRIPTION 

4.1  INTRODUCTION 

Computer  program  PLAN  is  developed  for  the  purpose  of  analyzing 
random  effects  in  a planar  array  of  thin-wire  dipoles.  It  is  written 
in  a general  form  with  several  options  included.  The  program  is 
especially  designed  for  statistical  analysis  but  it  can  also  be  used 
as  a general  program  for  analysis  of  deterministic  radiation  problems 
of  planar  arrays.  The  program  is  presently  dimensioned  for  a 10  by  10 
planar  array,  with  dipole  length  not  exceeding  A/2,  and  frequency  mode 
up  to  third  harmonic. 

The  program  first  reads  in  array  information.  Then  generalized 
impedance  and  admittance  matrices  are  calculated.  Next  is  the  procedure 
to  calculate  the  far  field  beam  patterns  for  a given  excitation  generated 
by  a random  number  generator.  This  procedure  is  repeated  many  times  to 
obtain  mean  and  standard  deviations  of  far  field  magnitude  and  gain. 

Two  disc  files  are  needed  as  auxiliary  storage. 


4.2  OPTIONS 

There  are  several  options  provided  in  the  program.  Each  of 
these  options  has  a corresponding  option  code  which  is  an  integral 


power  of  2.  The  options  are  specified  by  adding  up  all  the  code  numbers 
of  the  options  desired. 

(1)  1 - Generalized  Impedance  Matrix. 

(2)  2 - Compact  Admittance  Matrix. 

(3)  4 - No  Mutual  Effects. 

(4)  8 - Amplitude  Taper  and  Progressive  Phase. 

(5)  16  - Individual  Beam  Patterns. 

(6)  32  - Random  Phase  Errors 

Options  (1),  (2),  (5)  are  printout  options.  Option  (3)  represents 
a choice  in  the  method  of  solution  and  option  (4)  represents  a choice 
or  nominal  excitation.  (6)  is  a special  case  of  random  excitation  for 
which  the  phase  is  uniformly  distributed  between  -it,  it.  A brief 
description  of  each  of  these  options  follows: 

(1)  Generalized  ImpeddP';e  Matrix. 

The  generalized  impedance  matrix  [z]  for  the  array  may  be 
printed  out  by  selecting  this  option.  Because  of  the  block-Toeplitz 
property,  only  the  first  column  is  printed  out.  This  printout  option 
is  used  primarily  for  checking. 

(2)  Compact  Generalized  Admittance  Matrix. 

Recall  that  only  a portion  of  the  generalized  admittance 
matrix  [y]  = [z]  ^ is  computed,  namely  those  columns  corresponding 
to  the  non-zero  elements  of  the  excitation  matrix.  This  portion  of  the 
matrix  [y]  is  called  the  compact  generalized  admittance  ma  rix.  It 
may  be  printed  out  by  selecting  this  option.  The  compact  generalized 
admittance  matrix  is  k by  (M^M^M^)  where  the  latter  term  represents 


55 


the  number  of  unknown  expansion  functions  and  k represents  the  number 
of  dipoles  contained  in  this  rectangular  array  since  k columns, 
each  of  lengths  M2M3’  are  re<luire^  for  the  matrix  multiplication 
[y] [i] . Because  of  symmetry,  only  k/2  or  half  of  the  columns  are 
printed  out. 

Note  that  this  compact  matrix  contains  within  it  the  N-port  short 
circuit  admittance  parameters  of  the  antenna  system.  This  matrix  is  a 
kxk  square  matrix  within  the  compact  matrix  [22]. 

(3)  No  Mutual  Effects 

In  program  PLAN  mutual  effects  between  portions  of  each  di- 
pole and  between  different  dipoles  are  taken  fully  into  account.  If 
this  option  is  selected,  then,  in  addition,  the  program  is  run  neglecting 
mutuals  between  different  dipoles  but  taking  into  account  mutuals  within 
each  antenna.  This  additional  computation  is  extremely  rapid  as  com- 
pared to  the  full  analysis.  When  compared  with  the  complete  analysis, 
it  may  provide  some  insight  into  the  role  of  mutual  effects. 

(4)  Amplitude  Taper  and  Progressive  Phase 

If  this  option  is  selected,  then  the  user  may  specify  nominal 

* 

exponential  amplitude  taper  and  separate  prograssive  phase  shifts  in 

the  x and  z directions.  Thus,  the  beam  can  be  steered  to  an  arbitrary 
_ 

Note  that  amplitude  taper  function  can  readily  be  changed  by  replacing 
FUNX  and  FUNZ  with  any  desired  taper  function. 
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direction  in  (0,0)  space.  If  rhls  option  is  not  selected,  then  the 
excitation  is  automatically  made  uniform,  i.e.  all  non-zero  elements 
of  the  voltage  excitation  matrix  are  unity. 

A 

(5)  Individual  Beam  Patterns 

Note  that  the  mean  beam  pattern  magnitude  is  determined  by 
averaging  the  magnitudes  of  a number  of  separate  beam  patterns.  In 
some  cases  the  individual  beam  patterns  are  of  interest,  as  well  as 
the  mean  and  mean  minus  sigma  patterns.  If  this  option  is  selected, 
then  all  of  the  individual  beam  patterns  are  printed  out.  Some  examples 
are  shown  in  section  4.  This  option  is  of  interest  if  the  typical  beam 
pattern  differs  significantly  from  the  mean  beam  pattern.  This  is  the 
case,  for  instance,  if  the  random  phase  errors  are  very  large. 

(6)  Random  Phase  Errors 

This  option  selects  uniformly  distributed  random  numbers 
instead  of  normally  distributed  ones  for  the  phase  errors  of  the 
excitation.  This  option  is  of  interest  if  the  phase  error  of  the 
excitation  is  completely  random. 

4.3  MODELING 


(1)  Array  Orientation  - The  planar  (or  linear)  array  of  dipoles 
treated  in  this  program  lie  in  a plane  parallel  to  the  Z-axis.  For 
simplicity,  it  is  best  to  choose  an  xz  plane.  Each  dipole  is  identical 


and  z-directed.  Separate  uniform  spacing  is  required  in  each  dimension 
of  the  planar  array  (Fig.  8). 

For  simplicity  in  far  field  specification,  a planar  array  lying 


in  the  x-z  plane  with  dipoles  pointing  in  the  z-direction  is  recommended 
as  shown  in  Fig.  14 

(2)  Dipole  modeling  - All  dipoles  are  assumed  to  be  straight. 


Fig.  14.  A rectangular  planar  array  of  z-directed  dipoles  lying  in 
the  X-Z  plane. 


i 


58 


functions  used  per  dipole  is  determined  automatically  within  the  pro- 
gram. Three,  five,  and  seven  expansion  functions  are  used  for  L = A/2, 
A,  3A/2,  respectively. 


4.4  INPUT  DATA  SET-UP 


(1)  The  first  card  ontains  an  integer  number  which  is  the  summa- 
tion of  all  the  option  codes  selected.  See  section  4.2. 


Table  1.  CARD  FOR  OPTIONS 


Column 


Format 


Data  Description 


1-5 


Option  Code 


(2)  The  second  card  contains  information  on  the  frequency  mode. 
JF  = 1 means  fundamental  frequency;  JF  = 2 means  2nd  harmonic;  and 
JF  = 3 means  3rd  harmonic,  etc. 


Table  2.  CARD  FOR  FREQUENCY  MODE 


Column 


Format 


Data  Description 


1-5 


JF  - Frequency  mode 


The  JF  code  is  chosen  for  simplicity  in  treating  the  array  at  funda- 
mental and  harmonic  frequencies.  One  can,  of  course,  choose  any  fre- 
quency by  separately  choosing  antenna  length  L and  selecting  funda- 
mental frequency. 

(3)  The  third  card  specifies  the  array  data,  i.e.  the  number  of 
rows  and  columns  in  the  planar  array,  form  spacing  between  dipoles  in 


X Y Z-direction,  length  and  radius  of  each  dipole  (all  in  units  of 


wavelength) . 


Table  3.  CARD  FOR  ARRAY  DATA 


Column 

Format 

Data  Description 

1-5 

15 

M2  = if  of  columns 

6-10 

15 

M3  = if  of  rows 

11-20 

F10.5 

SX  - X-direction  spacing 

21-30 

F10.5 

SY  - Y-direction  spacing 

31-40 

F10. 5 

SZ  - Z-direction  spacing 

41-50 

F10.5 

HL  - Dipole  length 

51-60 

F10.5 

AK  - Dipole  radius 

(4)  The  fourth  card  contains  mean  and  standard  deviation  values 
for  both  amplitude  and  phase.  Also,  the  number  of  repetitions  K is 
specified.  If  mere  than  one  far  field  pattern  is  needed  for  one  exci- 
tation, set  the  pattern  code  JT  to  the  desired  number.  All  phase 
angles  are  given  in  degrees.  If  option  (6)  is  selected,  then  SD2  is 
ignored  by  the  program,  whether  or  not  a number  is  placed  there. 

Table  4.  CARD  FOR  VOLTAGE  EXCITATION  ERROR,  NUMBER 
OF  ITERATIONS  (K) , NUMBER  OF  FAR  FIELD 
BEAM  PATTERNS. 


Column 

Format 

Data 

Description 

1-10 

F10.5 

AMI  - 

Mean  value  of  amplitude 

11-20 

F10.5 

AM2  - 

Mean  value  of  phase 

21-30 

F10.5 

SD1  - 

STD  value  of  amplitude 

31-40 

F10.5 

SD2  - 

STD  value  of  phase 

41-45 

15 

IT  - 

if  of  repetition 

46-50 

15 

JT  - 

It  of  far  field  pattern 

(5)  If  option  of  amplitude  taper  and  progressive  phase  is 
selected,  then  the  desired  shift  angle  of  the  main  beam  in  both  di- 
rections must  be  specified  (in  degrees). 


(6)  This  card  specifies  how  the  far  field  beam  pattern  is  to  be 
selected.  All  far  field  points  are  assumed  to  be  equally  spaced  either 
a Long  a circular  arc  in  0 or  if  direction,  or  over  an  area  by  specifying 
the  starting  angle  of  CL  and  and  ending  angle  of  0^  and  <f>^.  The 

increment  angle  in  both  directions  is  specified  by  A0  and  A$  (Fig.  14). 
All  angles  are  in  degrees.  If  more  than  one  pattern  is  required,  then 
more  cards  specifying  far  field  points  are  needed. 
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Table  6.  CARD  FOR  FAR  FIELD  BEAM  PATTERNS 


Column 

Format 

Data 

Description 

1-5 

15 

NR  - 

Start  angle  of 

0 

6-10 

15 

NT  - 

End  angle  of  0 

11-15 

15 

NY  - 

Increment  of  0 

lt-20 

15 

NW  - 

Start  angle  of 

<f> 

21-25 

15 

NH  - 

End  angle  of  <p 

26-30 

15 

NJ  - 

Increment  of  <p 

' 
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5.  RESULTS 


i 

' 


Computer  program  PLAN  has  been  exercised  on  the  Honeywell  635 

computer  and  some  typical  results  are  shown  in  Figs.  15  -2Q. 

Fig.  15  shows  mean  and  m-o  beam  patterns  of  a 10  x 10  planar  array 

of  dipoles  (Refer  to  Fig.  1;  M2  = = 10)  for  various  values  of  0^,0^. 

The  nominal  excitation  is  uniform  and  the  nominal  pattern  (a  = a =0) 

a p 

is  shown  in  Fig.  15a.  0 represents  the  rms  departure  from  the  mean 

pattern.  It  is  assumed  that  the  nominal  voltage  excitation  is  1.0  volts. 
All  patterns  are  normalized  to  the  maximum  value  of  the  mean  pattern. 
Normalization  numbers  (N)  printed  on  each  plot  provide  a measure  of  the 
mean  field  intensity.  Note  that  the  mean  pattern  does  not  differ  signi- 
ficantly from  the  nominal.  The  range  of  variation  (0)  increases  with 

increasing  0,0.  The  effect  of  O is  more  pronounced, 
a’  p p 

Fig.  16  shows  the  corresponding  beam  patterns  for  the  second 

harmonic.  It  is  assumed  that  the  geometry  is  unchanged  from  that  of 

Fig.  15,  but  the  frequency  is  changed  from  to  2f^.  Since  the  spacing 

is  = 1.0  A at  second  harmonic,  grating  lobes  appear  at  <f>  = 0,  180°. 

The  nominal  pattern  (0  = 0 = 0)  is  shown  in  Fig.  16a.  Note  that  the 

a p 

mean  patterns  do  differ  more  significantly  from  the  nominal  than  was 
the  case  in  Fig.  15. 

Fig.  17  shows  the  corresponding  beam  patterns  for  the  third  harmonic. 
Note  that  the  grating  lobes  are  considerably  closer  to  the  main  beam. 
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The  nominal  pattern  (Q  = U = 0)  is  shown  in  Fig.  17a.  Note  that  the 

a p 

deviations  of  the  mean  from  the  nominal  pattern  are  significant  and  that 

the  range  of  variation  increases  with  0,0.  In  comparing  the  beam 

a p 

patterns  for  f^,  2f^,  3f^,  note  that  deviation  of  mean  from  the  nominal 
pattern  and  the  range  of  variation  (o)  both  increase  as  frequency  in- 
creases (for  given  0 , o ). 

a p 

In  addition  to  mean  and  (M-a)  patterns,  the  individual  patterns 
due  to  the  various  sample  excitations  are  also  of  interest.  The  computer 


program  includes  an  option  to  print  out  the  individual  beam  patterns 
from  which  t^he  mean  and  (M-o)  patterns  were  computed.  If  the  shape  of 


the  mean  and  (M-a)  patterns  are  relatively  close  to  that  of  the  nominal 
pattern,  then  this  option  may  not  be  necessary.  In  some  cases,  the  indi- 
vidual patterns  are  useful.  If  the  shapes  of  the  mean  and  M-o  patterns 
differ  significantly  from  the  nominal,  then  the  option  may  be  useful. 
Figs.  18  a,  b show  the  individual  beam  patterns  of  the  40  samples  used 

in  treating  a 5 x 5 planar  array.  In  this  case  (o  = 0.5,  o = 40°) 

a p 

the  sample  patterns  are  reasonably  close  to  the  mean  pattern  (Fig.  18c). 

Figs.  19  a,  b show  the  individual  beam  patterns  for  a greater  amplitude 

variation  ( o = 3.0,  0 = 40°).  Note  that  the  sidelobe  levels  differ 

a P 

greatly  from  the  mean  pattern,  but  that  the  main  beam  in  most  cases  is 
located  near  broadside.  Fig.  19c  shows  the  mean  beam  pattern. 

Fig.  20  a,b  show  the  individual  beam  patterns  for  a large  phase  variation 


(a  = 0.5,  0 = 180°).  In  this  case,  the  position  of  the  main  beam 

a p 

pattern  varies  considerably.  Fig.  20c  shows  the  mean  beam  pattern. 
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Note  that  ttu  individual  pattern  data  is  really  necessary  here.  The 
mean  beam  pattern  is  broad,  whereas  the  individual  beam  patterns  are 
narrow  but  with  positions  of  main  beam  varying. 

Fig.  21a  shows  the  individual  (sample)  patterns  of  an  eleven 
element  linear  array  for  the  extreme  case  of  completely  random  phase 
errors.  Note  that  the  beams  are  relatively  narrow  but  that  the  positions 
of  the  maxima  vary  considerably.  Fig.  21b  shows  the  mean  and  (m  - a) 
patterns.  Note  that  the  maximum  of  the  mean  beam  pattern  occurs  near 
broadside.  In  other  words,  there  is  a tendency,  despite  the  random  phase 
errors,  to  point  the  beam  near  broadside.  This  is  due  to  mutuals,  as 
shown  by  Fig.  20c,  which  indicates  the  mean  and  (m  - o)  patterns  of  the 
same  array  with  random  phase  excitation  for  the  special  case  of  no  mutuals. 
In  this  latter  case  the  direction  of  the  main  beam  appears  to  be  random. 


Computer  program  PLAN  includes  an  option  to  compute  the  beam 
patterns  without  taking  the  mutuals  into  account.  The  impedance  of 
an  individual  dipole  is  still  accurately  computed,  i.e.  mutuals  be- 
tween subsections  on  the  same  dipole  are  taken  into  account.  However, 
mutual  effects  between  dipoles  are  not  taken  into  account.  The  matrix 
inversion  then  becomes  very  simple.  Fig.  22a  and  b show  broadside 
patterns  of  a 5 x 5 planar  array  with  and  without  mutuals,  respectively. 
Figs.  22  c , d show  the  endfire  patterns  (array  scanned  exactly  to 
eidfire)  with  and  without  mutuals,  respectively.  Note  the  significant 
differences  in  the  sidelobe  levels  as  computed  with  and  without  mutuals. 

The  beam  patterns  of  a 9-element  linear  array  were  also  studied 


as  a function  of  o , o . Fig. 

a p 

array  as  a function  of  a , a . 

a p 


23  shows  the  sidelobe  level  of  the 
Note  that  sidelobe  level  increases 


monotonically  with  increasing  3,3  . 

a p 
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Figs  24-28  show  the  results  of  computation  of  some  of  the  gain 

quantities  G(j,  G^,  G^.  Fig.  24  shows  the  gain  of  a linear  array  of 

halfwave  dipoles  with  halfwave  spacing  (perpendicular  line  case)  as  a 

function  of  the  number  of  dipoles.  Computation  is  carried  out  with  and 

without  mutuals.  Note  that  the  gain  with  mutuals  is  always  larger. 

Figs.  25a  and  b show  the  ratio  G, /G.  as  a function  of  a , a for  an 

10  pa 

eleven  element  linear  array  of  halfwave  spacing  (perpendicular  line  case). 

rigs.  25a,  b represent  the  same  data  plotted  in  different  ways.  Note  that 

the  ratio  G^/G^  is  approximately  0.2  in  the  worst  case. 

Figs.  26a,  b show  the  ratio  G, /G_  as  a function  of  0 , a for  an 
n 10  pa 

eleven  clement  linear  array  of  halfwave  dipoles  with  0.6A  spacing  (parallel 
line  case).  Figs.  26a,  b represent  the  same  data  plotted  in  different  ways. 
The  data  of  figs.  25  and  26  differ  because  (1)  the  array  is  slightly  larger 
and  (2)  the  mutuals  are  different  in  the  parallel  line  case.  Note  that  the 
ratio  C^/G^  *s  ^ess  than  0.2  in  the  worst  case. 

Computations  of  G^/Gq  were  not  carried  out  for  a planar  array  because 
of  the  difficulty  in  locating  the  position  of  the  beam  maximum,  which  can 
occur  anywhere  in  space.  However,  using  the  array  factor  principle,  one 
may  deduce  approximately  the  corresponding  ratio  G^/ G^  for  an  11  x 11  element 
planar  array  by  multiplying  the  appropriate  factors  of  figs.  25  and  26. 

Thus,  the  worst  case  value  of  G^/G^  would  be  approximately  0.04  for  the 
11  x 11  planar  array.  This  represents  a gain  degradation  of  about  14  dB 


for  the  quantity  G. . 
• J. 
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Fig.  27a,  b show  the  ratio  G/G  as  a function  of  a , a for  an 

2 U pa 

eleven  element  1 inear  array  of  halfwave  dipoles  with  halfwave  spacing 
(perpendicular  line  case).  Figs  27  a,  b show  the  same  data  plotted  in 
different  ways.  Note  that  the  ratio  G^/G^  is  aPProximately  0.3  for  the 
worst  case.  This  would  imply  a ratio  of  0.09  or  a gain  degradation  for 
C'2  of  approximately  10.5  dB  for  an  11  x 11  element  planar  array. 

Figs.  28a,  b show  the  ratios  G^/G^  and  G2/Gg  as  3 ^unct^on  array 
size,  for  a linear  array  of  halfwave  dipoles  with  halfwave  spacing 
(perpendicular  line  case).  Computations  have  been  carried  out  both  with 
and  without  mutuals.  Note  that  the  gain  with  mutuals  is  larger  in  every 
case.  Note  that  G^/G^  decreases  to  very  small  values  in  the  random  case. 
The  corresponding  ratios  for  a planar  array  can  be  approximately  determined 
by  squaring  the  data  of  figs.  28  a,b.  For  a 3600-element  array  (60  by  60), 


gain  degradation  for  G^/Gq  is  in  excess  of  30  dB  and  gain  degradation  for 
G2/Gq  is  in  excess  of  20  dB  for  the  worst  case.  Thus,  gain  degradation  due 
to  phase  and  amplitude  errors  can  be  a significant  factor  in  determining 
harmonic  levels  in  the  far  field. 
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<f>  (DEGREES) 


ig.  15.  Mean  and  (M-<7)  broadside  beam  patterns  (in  the  principal 
H-plane)  of  a 10  x 10  planar  array  of  dipoles  at  funda- 
mental frequency  (L  = 0.5A,  S = 0.5A.  S = 1.0A 
a = 0.0005A)  x * z ’ 

(a)  a = a =0 
a p 
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Fig.  15.  Mean  and  (M-0)  broadside  beam  patterns  (in  the  principal 
H-plane)  of  a 10  x 10  planar  array  of  dipoles  at  funda- 
mental frequency  (L  = 0.5X,  S = 0.5X,  S = 1.0X, 
a = 0.0005X)  x z 

(b)  a = 0,  a = 

20° 

: 

a P 
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crp  = 40° 

N = 275.58 
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<t>  (DEGREES) 


Mean  and  (M-a)  broadside  beam  patterns  (in  the  principal 
H-plane)  of  a 10  x 10  planar  array  of  dipoles  at  funda- 
mental frequency  (L  = 0.5A,  S « 0.5X,  S = 1.0A 
a - 0.0005A)  x ’ z ’ 


Fig.  16.  Mean  and  (M-or)  broadside  beam  patterns  (in  the  principal 
H-plane)  of  a 10  x 10  planar  array  of  dipoles  at  second 
harmonic  (L  - l.QA,  S = 1.0A,  S - 2.0A,  a = O.OOiA) 

X Z 

(b)  a - 0,  ct  - 20° 
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Fig.  16.  Mean  and  (M-a)  broadside  beam  patterns  (in  the  principal 
H-plane)  of  a 10  x 10  planar  array  of  dipoles  at  second 
harmonic  (L  = 1.0A,  S = 1.0A,  s = 2.0A,  a = O.OOlA) 

X 2 

(c)  a *=  0,  CT  = 40° 
a p 


I 


0 30  60  90  120  150  180 

<f>  (DEGREES) 

Fig.  16.  Mean  and  (M-O)  broadside  beam  patterns  (in  the  principal 
H-plane)  of  a 10  x 10  planar  array  of  dipoles  at  second 
harmonic  (L  = l.QA,  Sx  = 1.0A,  S£  = 2.0A,  a = O.OOlA) 

(d)  ct  = 0.2,  a = 40° 
a p 
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Hg.  17.  Mean  and  (M-o)  broadside  beam  patterns  (in  the  principal 
H-plane)  of  a 10  x 10  planar  array  of  dipoles  at  third 
harmonic  (L  = 1.5A,  Sx  = 1.5A,  Sz  = 3.0A,  a = 0.0015A). 

(c)  a - 0,  CT  = 40° 
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Individual  and  mean  broadside  beam  patterns  (in  the 

principal  H-plane)  of  a 5 x 5 planar  array  of  dipoles 
at  fundamental  frequency  (L  =*  0.5A,  S “ 0.5A,  Sz  ■ 
1.0A,  a = 0.0005A,  o = 0.5,  a = 40°T 


(b)  Individual  (sample)  patterns 
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</>  (DEGREES) 

Fig.  18.  Individual  and  mean  broadside  beam  patterns  (in  the 

principal  H-plane)  of  a 5 x 5 planar  array  of  dipoles 

at  fundamental  frequency  (L  = 0.5X,  S = 0.5A,  S = 

1.0A,  a = 0.0005A,  o = 0.5,  a = 40°^f  z 

a 

(c)  Mean  and  M-a  patterns. 
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Individual  and  broadside  mean  beam  patterns  (in  the 

principal  U-plane)  of  a 5 x 5 planar  array  of  dipoles 

at  fundamental  frequency  CL  ■ 0.5X,  Sx  * 0.5X,  Sz  = 1 

a.  * 0.0005),  0 * 0.5,  <7  = 180°) 

3 P 

(a)  Individual  (sample)  patterns 
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Cb ) Individual  (sample)  patterns 
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Individual  and  broadside  mean  beam  patterns  (in  the 
principal  H-plane)  of  a 5 x 5 planar  array  of  dipoles 
at  fundamental  frequency  (L  = 0.5A,  S = 0.5A,  S = 1 
a = 0.0005A,  a = 0.5,  a = 180°)  x z 


(c)  Mean  and  (M-O)  patterns 


Fig.  21.  Individual  and  mean  broadside  patterns  (in  the 
principal  H-plane)  of  an  eleven  element  linear 
array  of  dipoles  at  fundamental  frequency 
(L  = 0.5X,  Sx  = 0.5X,  a = 0.005X,  a & - 0,  = “(random)) 

(a)  Individual  (sample)  patterns 


Individual  and  mean  broadside  patterns  (in  the 
principal  H-plane)  of  an  eleven  element  linear 
array  of  dipoles  at  fundamental  frequency 
(L  ■ 0.51,  S = 0.5A,  a = 0.005X,  a = 0,  a = “(random)) 

X cl  p 

(b)  Mean  and  (n-o)  patterns  with  mutuals 
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Mean  and  (M-a)  beam  patterns  (in  the  principal  H- 
of  a 5 x 5 planar  array  of  dipoles  at  fundamenta. 
frequency  (L  = 0.5A,  S = 0.5A,  S = 1.0A,  a = 0 
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Mean  and  (M-a)  beam  patterns  (in  the  principal  H-plane) 
of  a 5 x 5 planar  array  of  dipoles  at  fundamental 
frequency  (L  = 0.5X,  S = 0.5A,  S = 1.0A,  a = 0.0005X, 
0_  - 0.02,  a » 20°)  x z 


Broadside  patterns  - without  mutuals 


MUTUALS 


Gain  GQ  for  a uniformly-excited  linear  array  o 
dipoles  (L  - 0.5A,  Sx  = 0.5A,  a = 0.005A)  with 
without  mutuals  as  a function  of  the  number  of 
dipoles . 


Fig.  25.  Ratio  G^/Gq  for  an  eleven  element  linear  array  of 
dipoles  (L  * 0.5X,  Sx  = 0.5X,  a * 0.005X) 

(a)  as  a function  of  0 , a 
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Fig.  25.  Ratio  Gi/Gg  for  an  eleven  element  linear  array  of 
dipoles  (L  = 0.5A,  Sx  - 0.5A,  a - 0.005A) 

(b)  as  a function  of  a , a . 
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Fig-  27.  Ratio  G2/G0  for  an  eleven  element  linear  array  of 
dipoles  (L  = 0.5A,  Sx  = 0.5A,  a = 0.005A) 

(b)  as  a function  of  a . a . 
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Fig.  28.  Ratios  Gj/Gq  and  G2/G0  for  a linear  array  of  dipoles 
(L  = 0.5A,  Sx  = 0.5A,  a = 0.005A)  with  and  without 
mutuals  as  a function  of  the  number  of  dipoles 


(a)  OjGn 


6 APPLICATION  OF  THE  OTP  REGULATIONS 


6.1  THE  OTP  REGULATIONS 

Figs.  29,30  and  Table  1 show  some  of  the  current  regulations  [24]  of 
the  Office  of  Telecommunications  Policy  related  to  phased  arrays.  There 
are  several  salient  points  to  be  noted.  First,  at  the  band  edges  and  beyond 
^1  or  ^ ~ ^1^  » emissions  must  be  at  least  40  dB  down 

from  the  level  at  nominal  operating  frequency  Fq.  Second,  in  the  range 
f Fq  + Af 2 or  f £ ^0~^2  emissions  must  be  down  S dB  from  the  level  at 
fundamental  Fq.  S is  defined  in  Fig.  29.  Table  2 shows  typical  values  of 
Fq  and  P (maximum  power  spectral  level  in  dBm  /KHz)  for  which  S corres- 
ponds to  60  dB.  Third,  the  major  sidelobes  must  be  at  least  20  dB  down  and 
all  other  sidelobes  are  to  be  30  dB  down  from  the  beam  maximum. 

Revisions  of  the  OTP  Regulations  are  currently  being  considered.  Note 
that  there  is  no  explicit  restriction  on  the  maximum  level  at  Fq, 

however  S is  increased  in  direct  proportion  to  the  maximum  power  spectral 
level  P^.  In  other  words,  narrow-band,  high-power  systems  must  have  lower 
relative  levels  for  harmonic  and  spurious  emissions.  The  level  S may  be 
changed  somewhat  in  future  revised  regulations.  In  the  discussions  of  this 

chapter,  the  level  S is  assumed  to  be  about  60  dB,  corresponding  to  the 
cases  shown  in  Table  2.  Cases  falling  outside  this  range  will  be  discussed 
separately. 
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Fig.  30.  Weighting  Factor  k (OTP  Regulations) 


TABLE  1 


BANDWIDTH  DEFINITIONS  FOR  OTP  REGULATIONS 


„ 2kK  < 64 

B - - — (for  conventional  pulse  and  pulse  doppler  radars). 

= — + M ^ ■ + M (For  "modified"  pulse  radars  such  as  pulse 

compression,  chirp,  etc.). 


B = emission  bandwidth  in  MHz 


t = pulse  duration  in  microseconds 


t = pulse  rise  time 
r 

K = t/t 


BM  = bandwidth  due  to  pulse  modification  in  MHz 
k = weighting  factor  for  K (see  Fig.  30) 


B = 2Af x = Af2/5 


TABLE  2 


VALUES  OF  Fq, 

Fq  (MHz) 

100 

1000 

10,000 


\ 


Pt  CORRESPONDING  TO  S = 60  dB 

Pt  (dBm/KHz) 
0 

20 

40 
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How  does  one  demonstrate  fulfillment  of  the  requirements?  This 


question  has  not  been  completely  resolved  at  the  present  time.  Two 
possibilities  are  (1)  measurement  of  emission  levels  at  the  device  or 
at  antenna  terminals,  or  (2)  measurement  or  computation  of  emission  levels 
in  the  far  field  of  the  array.  Each  possibility  has  some  advantages  and 
disadvantages.  Technique  (1)  is  clearly  a simpler  measurement  (if  only 
a few  devices  or  terminal  pairs  are  tested),  whereas  technique  (2)  addresses 
a more  fundamental  quantity,  i.e.,  emission  levels  in  the  field.  Our 
results  show  that  in  general,  harmonic  levels  in  the  field  will  be 
significantly  lower  than  those  measured  at  the  device,  because  of  gain 
degradation  and  reflection  coefficient.  In  other  words,  the  array  has  a 
significant  effect  in  reducing  harmonic  levels.  Either  technique  should 
take  into  account  beam  scanning  and  operating  frequency  variation,  since 
Part  I [1]  of  the  final  report  has  shown  that  the  levels  of  fundamental 
and  harmonics  are  all  highly  dependent  on  load  conditions  which,  of  course, 
vary  with  frequency  and  scan  angle. 

Now  let  us  consider  some  of  the  points  noted  previously  concerning 
Fig. 29  . The  first  point  is  concerned  with  the  emissions  at  band  edges. 

The  emissions  at  band  edges  (F^  + Af^)  must  be  down  at  least  40  db  down 
from  the  fundamental.  Now  the  emission  characteristics  over  the  bandwidth 
B of  the  system  (F^  - Af^  - f - Fq  + Af^)  and  at  band  edges  depend  of  course 
on  the  details  of  the  particular  pulse  shape  used.  One  assumes  that  this 
requirement  can  be  met;  the  restrictions  which  it  places  on  the  radar  system 
will,  of  course,  depend  on  the  particular  system. 
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the  second  point  raised  concerns  the  levels  at  Fq  - Af 2 and  beyond, 
irtiich  must  be  down  at  least  S dB  from  the  fundamental.  This  requirement 
will  be  applicable  in  most  cases  to  the  2nd  and  3rd  harmonics  (2FQ  and 
3F0,  respectively).  For  instance,  the  bandwidth  of  the  bipolar  transistor  modules 
treated  is  less  than  12%.  Thus,  f^  would  be  less  than  0.06  FQ,  Af2  would 
be  less  than  0.6  F^  and  both  2nd  and  3rd  harmonics  would  lie  outside  the 
region  (FQ  - Af2  < f <_  F + Af  ) . In  these  cases,  then,  2nd  and  3rd 

harmonic  levels  must  be  at  least  60  db  below  fundamental  level.  This  require- 
ment is  discussed  in  sections  6.2  and  6.3. 

The  third  point  raised  concerns  the  sidelobe  level.  The  first  three 
sidelobes  must  be  at  least  20  db  down  and  all  others  at  least  30  db  down 
from  beam  pattern  maximum.  This  requirement  is  principally  a problem  of 
proper  array  design.  One  can  of  course  produce  sidelobes  which  are  suf- 
ficiently low  by  proper  amplitude  taper.  The  effect  of  random  variations 
in  amplitude  and  phase  can  be  predicted  using  the  computer  program  of  this 
report.  Then  the  taper  can  be  overdesigned  to  take  into  account  such 
random  variations.  The  variations  of  impedance  with  scan  angle  may  have  a 
significant  effect.  As  the  beam  is  scanned  and  the  load  presented  at  the 
output  of  the  solid  state  device  varies,  the  fundamental  output  will  vary 
as  pointed  out  in  Part  I.  These  effects  must  also  be  taken  into  account 
in  considering  the  desired  taper.  The  expected  variation  of  amplitude  and 
phase  could  be  used  as  a statistical  input  to  the  computer  program  of 
this  report.  In  other  words,  the  expected  variation  of  amplitude  and  phase 

could  be  used  to  determine  an  equivalent  0 and  O . Then  one  could  observe 

a p 

expected  variations  in  sidelobe  level  and  overdesign  accordingly. 
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In  summary,  the  sidelobe  level  requirement  is  considered  to  be  an 
array  design  problem  which  can  be  adequately  treated  by  the  usual  design 
techniques,  supplemented  by  the  computer  program  of  this  report.  In 
addition,  the  requirement  that  emissions  at  the  band  edges  be  40  db  down 
from  fundamental  is  considered  to  be  primarily  a pulse  design  problem. 

Thus,  our  attention  is  focussed  on  the  remaining  requirement  that 
emissions  outside  the  range  (F^  - Af^  - f - Fq  + Af2)  be  S db  down  from 
the  fundamental.  In  most  cases,  this  requirement  will  apply  to  the  2nd 
and  3rd  harmonics,  since  bandwidth  B will  in  most  cases  be  less  than 
0.20  Fq  and  F^  + Af2  will  thus  be  less  than  2Fq. 
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6.2  GUMMARV  AND  DISCUSSION  OF  EXPERIMENTAL  RESULTS 


Table  3 lists  approximate  levels  of  2nd  and  3rd  harmonics  of  the 
microwave  solid  state  devices  tested  with  matched  load.  In  parenthesis, 
the  range  of  levels  encountered  with  a VSWR  of  3.0  and  varying  angle  of 
reflection  coefficient  is  also  tabulated.  The  data  of  Table  3 is  extracted 
from  Table  6.1  of  Part  I of  Ref.  [1]. 

It  should  be  noted  that  several  additional  TRAPATT  devices  showed 
2nd  harmonic  levels  in  the  range  22-27  db  under  matched  conditions  (see 
Part  I).  The  reasons  for  the  superior  performance  of  the  type  3 circuit 
are  also  outlined  in  Part  I. 

For  the  moment,  we  will  assume  that  S = 60  dB  and  that  both  2nd  and 
3rd  harmonics  must  satisfy  this  requirement.  First,  consider  the  3rd 
harmonic  levels  listed  in  Table  3.  The  S-band  Bipolar  transistor  modules 
and  the  improved  (circuit  3 type)  TRAPATT  device  have  third  harmonic 
levels  which  are  lower  than  -60  db  relative  to  fundamental.  The  L-band 
modules  have  higher  levels  of  third  harmonic.  One  module  has  levels  at 
match  lower  than  -60  db  and  the  other  has  levels  at  match  of  -54  db.  The 
worst  case  is  -47  db  under  a VSWR  of  3.0.  These  third  harmonic  levels 
are,  of  course,  very  encouraging.  It  is  believed  that  very  little  effort 
has  been  made  to  reduce  2nd  or  3rd  harmonics  to  levels  as  low  as  -60  db 
and  it  therefore  seems  likely  that  with  improved  design  3rd  harmonic  levels 


can  be  reduced  to  less  than  -60  db  at  the  device  itself. 


i 

i 
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Now  consider  the  2nd  harmonic  levels  listed  in  Table  3.  Under 
conditions  of  match,  these  vary  from  -50  to  -54  db  for  the  L-band  modules, 

-37  to  -40  db  for  the  S-band  modules,  and  are  -51  db  for  the  improved 
(type  // 3)  TRAPATT  device.  The  L-band  modules  and  TRAPATT  device  are  thus 
within  better  than  10  db  of  -60  db  at  match  and  the  S-band  modules  are 
within  23  db  of  -60  db  at  match.  The  worst  cases  under  a VSWR  of  3.0  is 
-47  db  for  the  L-band  and  -34  db  for  the  S-band  modules.  Clearly,  these 
devices  require  some  improvement  in  regard  to  2nd  harmonic  levels.  It  may 
be  possible  to  significantly  lower  these  levels  by  design  procedures.  Some 
suggestions  are  made  in  Part  I for  appropriate  design  improvements.  Again, 
3ince  so  little  attention  has  been  given  to  lowering  harmonics  below  -60  db, 
it  seems  likely  that  some  improvement  can  be  obtained.  Since  the  worst 
case  for  L-band  modules  and  the  improved  TRAPATT  device  is  only  13  db  away 
from  the  -60  dh  level,  we  conclude  that  the  1-hand  modules  and  the  improved 
TRAPATT  device  (circuit  tvpe  *1)  are  strong  candidates  for  array  elements  in 
view  of  the  OTP  regulations.  Th->  vrorst  case  for  the  S-band  bipolar  modules 
is  27  dB  away  from  the  -60  db  level.  Olearlv,  some  significant  improvement 
is  reotiired  here.  Tt  should  be  noted,  however,  that  there  is  less  design 
experience  with  the  S-band  modules  than  with  the  L-band  nodules.  If  the  2nd 
harmonic  levels  of  the  S-hand  modules  cannot  be  improved  suf f icientlv,  it  mav 
he  necessary  to  add  a 2nd  harmonic  filter  to  each  module.  This  might  he  a 
reasonable  solution.  A microstrip  filter  should  not  greatlv  increase  cost  and 
size. 
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In  conclusion,  the  2nd  and  3rd  harmonic  levels  are  encouragingly  low 
for  most  devices.  The  -60  db  levels  are  not  satisfied  in  all  cases  at  the 
device  itself.  However,  many  of  the  levels  are  close  to  the  -60  db  levels 
required  and  it  seems  likely  that  modest  improvements  in  design  should  permit 
most  of  these  devices  to  satisfy  OTP  regulations  at  the  array  terminals. 

The  S-band  bipolar  modules  require  significant  improvement  in  2nd  harmonic 
levels  and  an  additional  filter  may  possibly  be  required.  Thus,  we  conclude 
that  of  the  solid  state  devices  tested,  the  L-band  bipolar  modules  and  the 
improved  TRAPATT  devices  are  strong  candidates  for  array  elements  and  that 
the  S-band  bipolar  modules  are  also  viable  candidates. 

For  those  radar  systems  for  which  S <_  60  dB,  then,  one  expects  that  with 
some  improvement  in  design,  the  OTP  regulations  can  be  satisfied  at  the  antenna 
terminals.  In  addition,  a significant  gain  degradation  is  expected  at  harmonic 
frequencies,  as  noted  in  chapter  5. 

A point  should  be  made  here  concerning  those  radar  systems  for  which 
S > 60  dB.  For  such  systems,  the  design  problems  and  especially  the 

measurement  problems  may  be  quite  severe.  However,  the  additional  effect 
of  the  array  in  reducing  harmonic  levels  may  be  predictable  by  using  the 
program  of  this  report,  and  in  some  cases  this  effect  may  be  sufficient  to 
reduce  the  emission  in  the  field  to  the  required  levels.  The  effect  of 
the  array  is  discussed  in  Section  6,3, 
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6.3  EFFECTS  OF  THE  ARRAY  ON  FAR  FIELD  LEVELS  OF  HARMONICS 


! 

1 


The  result  of  computation  of  Chapter  5 have  shown  that  the  array  has 
a significant  effect  upon  harmonic  levels  in  the  field.  In  particular,  if 
phase  and  amplitude  are  highly  random  at  harmonics,  then  there  will  be  a 
significant  degradation  of  gain  as  indicated  by  figs.  24-28. 

In  addition  to  the  gain  degradation,  an  array  of  halfwave  dipoles  will 

usually  be  well  matched  at  fundamental  F^,  highly  mismatched  at  2nd  harmonic 

2F0>  and  fairly  well  matched  at  3rd  harmonic  3F^.  Fig.  31  shows  the 

2 

reflection  factor  1 - |r|  (where  T is  the  reflection  coefficient  at  antenna 
terminals),  for  a dipole  of  length  L and  diameter  D,  fed  with  a transmission 
line  of  characteristic  impedance  of  50  ohms.  The  data  were  computed  by 
the  method  of  moments.  Note  that  there  is  a significant  reflection  loss, 
of  the  order  of  10  dB  or  greater  in  some  cases,  at  2nd  harmonic  (L  = 1.0A). 
The  eftect  of  this  reflection  factor  is  complex.  The  load  presented  to  the 
solid  state  device  at  F^  is  different  from  that  at  2F^  and  one  can  determine 
the  corresponding  harmonic  levels  only  by  test,  in  general.  In  some  cases, 
however,  experiments  have  shown  [1]  that  the  levels  generated  depend  pri- 
marily on  the  load  presented  at  fundamental  frequency.  In  such  cases,  one 
can  then  use  t.ie  reflection  factor  at  2nd  harmonic  directly. 

Thus,  there  are  three  principal  effects  of  the  array  on  harmonic  levels 
The  first  is  a degradation  of  gain  due  to  phase  and  amplitude  errors  at 
harmonic  frequencies,  the  second  is  an  increased  electrical  array  size  at 


2i 
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Fig.  31.  Reflection  Loss  for  a centerfed  dipole  of  length  L, 
diameter  D,  fed  by  a 50  ohm  transmission  line,  as  a 
function  of  antenna  length. 


7.  CONCLUSIONS 


R 1 


The  application  of  the  OTP  Regulations  to  Solid  State  Arrays 
has  been  considered.  In  part  I [1]  of  this  final  report,  typical 
emission  levels  of  some  solid  state  devices  are  reported.  Using  this 
data,  as  summarized  in  table  3 of  part  II,  it  is  concluded  that  the 
devices  tested  (Bipolar  transistors  and  TRAPATT  devices)  are  reason- 
able candidates  for  array  elements,  assuming  that  the  level  S is 
approximately  60  dB.  For  S levels  in  excess  of  60dB,  significant 
development  and  experimental  problems  are  anticipated  if  measurements 
are  to  be  made  at  the  antenna  terminals.  However,  the  array  itself 
is  significant  help  here  because  of  gain  degradation  and  reflection 
loss  at  harmonic  frequencies.  A computer  program  has  been  prepared 
to  analyze  random  effects  in  planar  arrays  of  dipoles.  Typical  results 
indicate  that  the  array  has  a significant  effect  in  reducing  the  ratio 
of  harmonic  to  fundamental  levels.  If  this  factor  is  taken  into  account 
the  solid  state  devices  become  even  stronger  candidates  for  array  ele- 
ments in  view  of  the  OTP  regulations. 
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Appendix  A 

The  Matrix  Inversion  Routine  for  the  Impedance 
Matrix  of  a Planar  Array 

The  algorithm  used  in  this  program  to  obtain  the  inverse  of  a 
block-Toeplitz  matrix  is  based  on  D.H.  Sinnott's  implementation  [19] 
and  can  be  summarized  as  follows : 

For  a given  block-Toeplitz  matrix: 


-zi 


Z1 -ZN 


"N-l 


[Z] 


N-l  ...Z 


0 - 


where  each  submatrix  Z.  is  of  order  n. 

-~i 

be  partitioned  in  the  same  fashion. 


Then  inversion  matrix  [Y]  can 


where  are  square  matrix  of  order  n. 

The  computer  storage  requirement  for  this  algorithm  is: 

Z((N  + 1)  x n2) 
ip(2N  x n2) 

A(2n2) 

Y(^(N  + l)2  x n2) 

Since  matrices  ip^  are  intermediate  data,  and  are  calculated  itera- 
tively in  Eq.  (1)  to  (4),  i.e.  \p^2\  i^2) » 

they  can  be  stored  alternately  in  two  disc  files  A and  B instead  of  in 


)-»■  disc  B 


This  saves  the  storage  requirement  for  the  tp  array.  In  addition,  since 
submatrices  Y are  only  a function  of  A and  ip  in  Eq.  (5)  and  (6) , it 
can  be  restored  into  the  Z array  if  [Z]  and  [Y]  are  not  subsequently 
required.  Due  to  the  effects  of  the  zero's  in  the  voltage  excitation 
matrix  as  explained  in  section  2.6,  only  part  of  the  [Y]  matrix  is 
computed,  i.e.  a condensed  [Y]  matrix.  Hence,  the  size  of  the  array 
[Y ] is  almost  compatible  with  the  Z array,  which  saves  the  storage 
requirement  for  [Y]  array. 

The  resultant  program  needs  computer  storage: 

Z( (N+l)  x n2) 

A(2n2),  ST(n2),  PSl(n2),  PS2(n2) 

where  ST,  PS1,  PS2  arrays  are  used  as  temporary  storage.  The  computer 
storage  requirement  is  cut  down  to  only  a small  fraction  of  the  original 
requirement. 


In  Eq.  (5)-(6),  submatrices  Y are  calculated  subsequently  from 

Y , , , A,  and  \p.  Since  there  are  only  a few  fixed  columns  of  each 

-r-l,s-l  - Z 1 

submatrix  Y required  with  respect  to  non-zero  element  in  [V],  once 
rs 

the  fixed  columns  Y and  Y.  submatrices  are  obtained  for  0 < r < N, 
rO  Or  — — 
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the  fixed  columns  of  all  other  submatrices  Y can  be  determined  without 

-rs 

calculating  the  other  elements  of  the  submatrix. 

1 2 3 

The  number  of  complex  multiplications  is  approximated  as  y(N+l)  (n) 
for  Eq.  (1)  - (A)  and  y (JH-l)^(n)^  for  Eq.  C5)-(7).  When  taking  ad- 
vantage of  the  condensed  [Y]  matrix,  the  number  of  complex  multiplications 

1 3 

for  Eq.  (5)-(7)  is  about  y(N+l)n  , which  is  relatively  small  compared 
1 2 3 

with  y (N+l)  n , hence  it  can  be  neglected.  The  resultant  program  needs 
about  one-half  the  complex  multiplication  required  by  the  original  one. 

A typical  example  of  comparison  between  three  matrix  inversion  routines 
is  shown  below. 

For  a block-Toeplitz  matrix  of  size 

(N+l)xn  = 1400  - order  of  the  entire  matrix 

N+l  = 20  - order  of  the  matrix  in  terms  of  blocks 

n = 70  - order  of  submatrix 

Computer  memory  required  is  listed  below  (in  units  of  million  words) 


Gauss-Jordan 


Sinnott 


Improved 


The  number  of  complex  multiplications  (in  units  of  million)  ere: 


Gauss-Jordan 


Improved 
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COMPUTER  PROGRAM  - FLAN 
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- BY  PETER  HSI 

- NOV.  10#  L975 

* ANALYSIS  OF  THE  FIELD  PATTERN  AND  GAIN  DEGhEGATI ON 

* OF  A GIVEN  PLANAR  ARRAY  DUE  TO  RANDOM  NOISE  EXIST 

* ON  EACH  DIPOLE 

•DEFINITION-  M 1 = #0  F EXPANSION  FUNCTIONS  ON  EACH  DIPOLE 
M2=#0F  COLUMNS  IN  PLANAR  ARRAY 
M3*#0F  ROWS  IN  PLANAR  ARRAY 
HL=DI POLE  LENGTHC IN  WAVE  LENGTH) 

AK=DI POLE  RADI USC IN  WAVE  LENGTH) 

SX=X- DI KECT ION  SPACING  BETWEEN  DIPOLES 
SY *Y-  * " 

SZ  = Z-  “ " 

JF=FRFQUENCY  MODE#  1 - FUNDAMENTAL 

2-SECOND  HARMONIC 
ETC. 

NW»MAX. (M2#M3) 

NP=MI  N . ( M2#  M3  ) *M  1 

NNN= TOTAL  # OF  FAR  FI  ELF  POINTS 

•STORAGE  REQUI REMENT- 

Z(NW*NP**2)  OR  Z( 0. 5*M1*CM2*M3)**2) 

CUR#  CUR2#  CC# DD(M1*M2*M3) 

HX,HY#HZ, V, VP# (M2* M3) 

TE#  TEE#  PH#  PHH,  T0T1#T0T2»T0T3#T0T4#  (NNN  ) 

A. A#  BP  ST1#ST2,  ST3#  STA,  < NNN  ) 

Z2(M 1**2) 

TWO  DISC  FILFS  FOR  TEMPARARY  STORAGE 


COMPLFX  Z ( 5000) » Z2( 50) 

COMPLEX  ZZ 

COMPLEX  CUR(300)#CUR2(300),V( 100)#VP( 100) 

DIMENSION  HX( 100), HY( 100)»HZ( 100) 

DIMENSION  TOTH  1S0)#T0T2(  1S0)»T0T3(  180)#T0T4(  180)# 
5TE(  180)#PH( 180)#TEE(  180)#PHHC 18  0)# ST  1C  180),ST2C 180) 
DIMENSION  ST3(  18  0) # STAC  180) # G 1 C 1 00) # G2C 1 00) 

COMMON  AAC 180)#PP( 180) 

LOGICAL  I JF(  1 0) # Y 

DIMENSION  Y2<  180)  # YMINC  3)  , YMAX  ( 3 ) # YCHARC  3)  , Dl'MY{  180 

COMPLr X CCC300)#rT(300),FE 

DATA  (YCHARC I ) # I a 1 # 3 ) / 1H* # 1 Hf # 1H#/ 

KI  =07 

CALL  FLCF0F(5#Y) 

CALL  OFTION(IOP) 

CALL  Ai.x.AY  ( MX  • HY  • HZ  . HI..  AK.  D?  .M  1 . MP.M3) 
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M23=M2*M3 

CALL  IMPDCM1,M2,M3,DZ, AK»  HX,  HY,  HZ,  Z > 

NW=M3 

IFCM3.LT.M2)  NV=M2 
NB=M2*M3/NW 
NP=M l*NB 
L = 0 

K=  C M 1/  2) *NP 
DO  82  I = 1,NW 
DO  83  J=  1 * NP 
L=L*  1 
K=K*1 

CCCL)=ZCK) 

CONTINUE 
K=K*NP*  CNP-  1 ) 

I FC IOPC 1 ) ) GO  TO  111 
KKK=NP»NP*NW 
PRINT  61 

FORMATC/// IX, 'GENERALIZED  IMPEDANCE  MATRIX',//) 
PRINT  60, C Z C I >, 1= 1, KKK) 

FORMAT C IX, 10F1 1.6) 

CONTINUE 

CALL  LINVRTCZ,NW,NP,NB) 

NH=  CNW+ 1 )/2 

MP=NP*NB 

MPN=MP»NW 

NPT=NH*NB*NV*NP 

IFCI0PC2))  GO  TO  100 

PRINT  91 

FORMATC/// IX,  'COMPACT  ADMITTANCE  MATRIX',//) 

PRINT  60,  (Z( I ), Is 1,NPT) 

CONTINUE 
REWIND  KI 
DO  90  LL= 1 , NH 
DO  90  J = 1,NB 
DO  90  1=1, NW 
DO  90  K=  1 , NP 

L=  Cl  - 1 ) “MP*  C LL-  1 > *MPN+  C J-  1 ) *NP*K 
ZZ=Z(L) 

VRITECKI)  ZZ 
CONTINUE 

IFCIOPC3))  GO  TO  300 
CALL  N OMUL CM  1, DZ,  AX,  Z2) 

NN2=M 1*M  1 
PRINT  71 

FORMAT C /// IX,  'SELF- I MPFDAVCE  C.\0  MUTUAL  EFFECT)'//) 
PRINT  60,  C Z2C I ) , I = 1 , NN2) 

CALL  LIN  EOC  Z2, M 1 ) 

PRINT  81 

FORMAT C /// IX,  'SELF- ADMI TTANCF  CNO  MUTUAL  FFFFCT)  '// ) 
PRINT  60,  C.Z2CI  ),I=  1,NN2) 


"”"T'  ’ IHf  ■«» 


r I 


dc  23  i = i > r:ri2 

33  Z2C  I ) =7.  ?A  I >*0.  0 1 

300  emdfilz  :<i 

L = 0 

u I23=:n*ti2*ti3 
K=  0 

DC  9 3 1=  1/NV 
DC  93  J=l/iJP 
L=  L+  1 
K = K+  1 

93  DDC  L) =7 (K) 

9 2 K=1C+MP*CNE-  1 ) 

ee=cmplxco . / o . ) 

DO  9 4 I = Ml  23 
EE=  EZ+CC ( I ) *DD  C I ) 

9/1  cc:jti  jt'e 

FEINT  9 6/ EE 
TTP I TEC  2/  9 6 ) EE 

FO°:iAT(  ///lYz  'CHECK  FOP  MATP.IY  INPEPSIOM  '/2F12.6//) 
TP=  6*2831853 
PE'/  I HD  K I 
DO  66  1=1/NPT 
r EAD  ( IC I ) 7.7. 

Z < I > = 7.7*0  . 0 I 

PEAD  40/AUUAM2./  FD 1 / 5D2/  17/  JT 
I F(  v ) CALL  EYI7 

FGP;iATC///l  Y.*'  it  CF  I TEFATIOHS  = % I 10*//) 

T'TC  IT  EC  2/  19)  AMU  A!  12*  SD1/  SD2 
P P I 'IT  7 0 7 

fo  r\:xr(.//j  if*  'raiFcru  e:'CI'"atic:i  applied  -**//> 

P P I M T 1 9*  Ail  1 * AU2/  FD 1 * SD2 
?D2=3D24>T°/360. 

F CP  fiATCAFlO. 5/215) 

PPI1JT  823/  IT 
DC  250  I=l/M23 
’'PC  I )=C;!rLv(  1 . 0/  0 • 0) 

I FC I 0~C/i)  ) GO  TC  30  1 
CALL  7 A”  EP  C * rP  / U 1 / il  2 / I-i  3 / BY  / HY  / H Z ' 

COM  T I i!TTE 
DC  2C1  J=1/JT 
l'/  - 1 35  79 

P ZAC  C C / Mr / NT / MY zIJK/HS/MJ 
FCT,MA'r<  615) 

A ’ ' 3 1-3  • 

A’/G2=0. 

DO  5 C 1=1/130 
T 0'r .1  C I ) = 0 . 0 
T CT2C  I ) = 0 • 0 
7 CT3C  I ) = 0 • 0 
'**  CT'A1  ) = 0-  0 
~EZ<  I)  = 3.  9 

i) -c • c 
7 e c i ) = c . 0 
I ) . c 
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DC  1 000  K=  1 / I T 

CALL  <raLT<'^^P>M£*M3>  IX*  AM I*  AM  2*  SD1*.SD2> 

CALL  H T 'LT  ( £ * V,  C nr * 11 1 * H 2 * ! 1 3 > 

CALL  FAr < CUP*  nr * MT  * HY*  MK* H S*  N J * HX* HY * H Z * NNN * DZ  * 

1PH*TE*T0T1  *T0T2*I123*!ll*  V*  AVG  1 ) 

I F(I  cri  5)  ) GO  TO  101 

PTUMT  13*K 

N A=Mr>+  1 

N E=NY+ 1 

NC  = .MK+1 

MD=MS+1 

LL=  0 

El  G=  0. 

DO  11  LY=MC*  ND*NJ 
DC  11  LK=NA*  HE*  NY 

LL=LL+ 1 

ith«lx-i 

IPC=LK- 1 

PT’l:r  15*1  TH*  I PC  * AA. ( LL)  * DEC  LL ) 

I FCAACLL) .GT.EIG)  BIG=AA(LL> 

11  COM  TXfJUE 

DC  333  1=D3 

Y MAY  (I ) = 1 . 0 

333  Y .11 ']  ( I)  = 0 • 0 
DO  334  1=1* NNN 

334  AA  ( D =AAC  I ) /BI  G 
P^IM  T 18*EIG 

CALL  ge^lotc  o*.\’Bn*  0*  0 • * 3»*  YCKAP*  YMAX*YiiIU*  AA>  DPMY* 

£ DTJ  MY  * DtL-lY  * DUMY  * DUMY  > 

101  CC'l  T I MV E 

I FCI  C-1  (31)  GO  TO  1000 

CALL  MVLT3(T2*”*  C"°£*H  I*:i2*:13> 

CALL  F AF  C CUT  2 * ’ T * ; T * ; J Y * 1 1 1 * M S > H J * 1 D * HY  * H Z * ! 1 M M * 

S D Z * ^ ; i:  i * T E E*  T CT  3 * T C ' T 4 * i 1 2 3 * M Ij  V * A ' VG  2 ) 

1 0 00  COM  T I MUE 
PFIMT  166 

166  FO  FMAT( ///IX*  ' FA?  FIELD  PATTEPN  ( E- FI  ELD)  '*  //) 

CALL  CTDF(TOT  I * T CT2*  PH* TE*  ME*  MT*  MY* HK* M S* MJ*  STl*  ST2*  IT*A"G1  > 
IF aP(o>)  GC  TC  177 
PEI MT  167 

167  FC?  1ATC ///IX*  * FAr  FIELD  PATTEPM  ( E-FI  ELD)  '* 

S'  **»  tie  MTH-T'AL  EFFECT  TAKIMG  IMTC  ACCC"’JT  '//> 

CALL  STDV(TCT3*?CT4*  "HH* TEL*  !1°* MT*  : Jv,  ,'JK*  MS*  M«J*  ST 3/ 

SST4* IT*A’'G2> 

17  7 N C=  0 
1 7C  CCMTIWUE 

DO  110  1=1*  3 
YflA:;<I)  = -l  .E-36 
110  Ylili-KI  > = 1 .E36 
DC-  120  1 = 1*  MMM 
Y2(  I >=TOT  t ( I )-STl  < I ) 

Y MAY  < 2)  =Ai  1AX  1 (Y.1AX  ( 2 > * TOT  1(D) 
v.;i;;(2)  =AMI  M 1 ( Y.I IM  < 2)  * TOT  1 ( I > ) 

120  C C.l TIMED 
20  DC  170  1 = 1/  M MM 

Y2( I ) =Y2( I ) /YMAXC2) 

17  0 TCT  1 ( I ) = TAT  1 ( I)  /Y.-lAY  ( 2 ) 
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A 





- ■ 


WRITE  (6*  19)  AM1.AM2*  SD1*SD2 
WRITE  !6*  16) 

DO  160  1 = 1*  2 
YMAr<I)  = l .0 

160  y:ii:jci)  = o .0 

CALL  GEPLOT!  2*NNN*  0/  0w3.  * YCHAP* YMAX* YMIM*  TOT  1 <Y2< 

& DUMY  * DU  MY*  DU11Y  * DUMY ) 

IFCM  C.EO.  1)  GO  TO  201 
I F(I0p(3)  > GO  TO  20  1 
DO  7 7 I=1*NNN 
TOTl!  I)=T0T3!I > 

ST  1 ! I ) =ST3( I ) 

7 7 CONTINUE 
N 0=  1 

GO  TO  178 

10  FORMAT!  2X  * 13F6.3//) 

12  FORMAT ( 61  12) 

15  FORMAT!  1Y*2I 10* IR2E30.6) 

13  FORMAT  (// IX*  'INDIVIDUAL  PATTERN  **'*110*  //*  8X* 

1 'THE'*  7X * 'PHI  '*2 OX*  'E-FIELD'*  20X*  'CAIN  '*  //,) 

14  FORMAT C 2X  1P3E14.5) 

16  FOR  MAT  ! 1 H 1 > 

9 FORM  AT  (6F 12. 6) 

18  F CPMAT ! IK  1 //  22HNCPMALI ZATI ON  FACTOR  = 1PE14.6) 

19  FO  PMATC //12HMEAN  VOLTAGE  9X  1H=  1 PE  1 4 . 6/ l 0HMEAN  PHASE 

1 1 IX  1 H=  1PE14.6/1 6H VOLTAGE  VARIANCE  5X  1H=1PE14.6/ 

2 14HFHASE  VARIANCE  ?X  1H=  1PE14.6) 

201  CONTINUE 

30  GO  TO  7 
END 


SUBROUTINE  ARRAY! HX»  HY » HZ  * HL»  AK*  DZ * e.  i » M2* M3 ) 

C READ  IN  ARRAY  DATA 

DIMENSION  HX ( 1 ) * HY ( 1 ) * HZ ( 1 ) 

TP=6. 2631853 

READ  10,M2»M3,SX,SY,SZ,HL,AK 
READ  1 O#  JF * M 1 
10  FORMAT C 21 5*  5F1 0. 5) 

M23=M2»M3 
PRINT  30, M2. M3 
WRITE!  2,  30)  M2, M3 

30  FORMAT!/// IX,  ' ******  PLANAR  ARRAY  PROGRAM  *****',///, 
$1X,  ' PLANAR  ARRAY  OF ',15,'  BY  ',12,'  DIPOLES',//) 

LL*0 

DO  20  1*1. M3 
DO  20  J* 1 , M2 
LL*LL* 1 

HX!LL>»! J- 1)*SX 
HY!LL>»! J- 1)*SY 
HZ!LL)*!I-l)#SZ 
20  CONTINUE 
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PRINT  112 

FORMAT C 3X»  'ELFMFNT  '»  1 5X,  'LENGTH ',  1 2X,  'X-COORDIN  ATE ',  8X, 
1 'Y-C00RDINATE',8X,  'Z-COO  RDI  NATE ',  8X,  'RADIUS',///) 

PRINT  1 1 1,  Cl »HL»HXC I ) ,HYC I ) , HZ C I ) , AK, I = 1 , M23 ) 

FORMAT C 5X, I 3,  8X, 5E20. 7) 

DO  A3  1 = 1#  M23 

HXC I ) »HX ( I )* TP* FLO AT C JF) 

HYC I ) =HYC I )* TP* FLO AT C JF) 

HZ  Cl )*HZC I )*TP*FLOATC JF) 

CONTINUE 

HL=HL*TP*FLOATC JF) 

AK=AK*TP»FLOATC JF) 

DZ*HL/FLOATCMH-l> 

RETURN 

END 


SUBROUTINE  OPTIONCIOP) 
DETFRMINF  OPTIONS  BEING  CHOSEN 
LOGICAL  IOP(l) 

REAP  10, I OPT 
FORMAT C I 5) 

DO  20  N= 1 , 10 
K*1  I-N 

INDEX* 2**  C K- 1 ) 

IOPC  K)  = . TRUE* 

IFCIOFT.LT. INDEX)  GO  TO  20 
IOPT=!OPT-IN  DEX 
I0PCK)=. FALSE. 

CONTINUE . 

RETURN 

END 


SUBROUTINE  IMPDC M 1 , M2,  M3,  EZ, AK,  HX,  HY,  HZ,  Z ) 

GENEhATE  IMPEDANCE  MATRIX  DY  DLOCK-TOEPLI TZ  PROPEhTY 
- RETURN  WITH  SUPMATRIXES  Z 0, Z 1 , Z2, . . . ETC. 

DIMENSION  HX  Cl), HY  Cl), HZ  Cl) 

COMPLEX  AC  140), ZC 1 ),ZMNG 
11=0 

IFCM3.LT. M2)  GO  TO  60 

M 1 2*M 1 *M2 

DO  30  M=1,M3 

MM= 1 ♦ C M- 1 )*M2 

DO  31  L=1,M2 

LL=CL- 1 > *M 1 

ZM-HZCL) 

ZN-HZCMM) 

SD1«HXCL)-HXCMM) 

SD2*HY  C L ) -HYCMM ) 

SD3  = HZ  C L ) - HZ  C MM ) 


SD«SORTC  SD1*SD1  + SD2*SD2+SD3*SD3> 
IF(SD.EG.O.)  SD=  AK 
DO  21  I»1,M1 

ACLL+I >»ZMNGCZM,ZN, DZ,  DZ,  SD) 

ACLL+I  >«A<LL*I  )*0.01 

ZM*ZM*DZ 

CONTINUE 

IF(M.EG.l)  GO  TO  22 
DO  23  L«1,M2 
LL=CL-1)*CM1-I> 

KK=MM+L-1 
ZM-HZC l> 

ZN«HZ(KK)*DZ 
SDl-HXC l)-HX(KK) 

SD2-HYC  1 )-HY\'KK) 

SD3*HZ ( 1 >-HZCKK> 

SP=S0RTCSD1*SDUSD2*SD2+SD3*SD3> 

DO  24  1=2, Ml 

AC  M 1 2+LL* l - 1 ) =ZMNGC ZM, ZN , DZ,DZ,SD> 

A(M18*LL*1  - 1 )*A(M't2*-LL+l  - 1 >*0»  01 

ZN»ZN^DZ 

CONTINUE 

CONTINUE 

MMM-M 

MM3-M3 

CALL  FLKCM 1 , M2, MM3, MMM, A,  Z> 

CONTINUE 

GO  TO  70 

DO  40  M»  1 , M2 

LL-0 

DO  41  L* 1 , M3 
L1=M*CL- I)*M2 
ZM=HZ( LI) 

ZN»HZC  1) 

SD!=HX<L1 >-HXC 1) 

SD2*HY ( L 1 > -HY  C 1) 

SD3=HZ(L1)-HZ( 1) 

SD=SQHT(SDl*SDl«-SD2*SD2-*-SD3*SD3> 
IF<SD.EQ.O.)  SD»AK 
DO  42  1=1, Ml 
LL-LL+ 1 

A(LL)=ZMNG<ZM,ZN,DZ,DZ,SD>*0.01 

ZM-ZM+DZ 

CONTINUE 

rFCM3.EC. 1 ) GO  TO  90 
DO  50  L3?, M3 
L1=M«.<l-  1)*M2 
ZM=HZ( 1 ) ♦DZ 
ZN=HZCL1> 

SD1*HXCL1 >-HXC 1) 

SD2-HYCL1 >-HYf 1) 

SD3=HZCL1)-HZ(  1) 

SD=SQRTCSD1*SDUSC2*SD2^SD3*SL3> 
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ihMihmi 


_ 


IFCSD.EQ.O.)  SD«AK 
DO  51  I =2,  Ml 
LL-LL*  1 

ACLL)=ZMNGCZM,ZN, DZ, DZ, SD> 
AC LL ) * A(  LL) • 0»  0 l 
51  ZM-'ZM*DZ 
50  CONTINUE 
90  MMM-M 
MM2-M3 

CALL  SLK2CM I * MM2,  MMM,  A,  Z > 
AO  CONTINUE 
7 0 RETURN 

END 


SUBROUTI  NE  FLKCM  1 , M2,  M3.  M,  A,  R) 

C GENERATE  A TOEFLITZ  MAI  hi  X PY  1-ST  HOW  $ COLUMN 

COMPLEX  AC  1 >,EC  1 > 

M 1 2*M 1 *M2 
M 1 23*M 1 2*M 1 
M123A»M123*M2 
DO  »00  N = 1 , M 2 
DO  100  1*1,  M2 
DO  100  J = 1 , M 1 
DO  100  L= 1 » M 1 

K1=C AP5CL-J0+ 1 ) *M 1 * APS  C I - N ) 

IFCM.EC.  1 ) GO  TO  1 5 

IFCL.LT. J)  K1=K1  + M12-CAPSC I -NO*  1) 

1 5 K»L+M  1*C  I -1)*M12*CJ-  1)+-M123*CN  - 1 )*M1234*CM-  1 ) 

PCK)=ACK1 ) 

100  CONTINUE- 

DO  1 10  K=  1,M12 
I *K+  CM- 1)*M123A 
110  CONTINUE 
12  FORMAT C 2X, 9F6. 0) 

RETURN 

END 


FUNCTION  ZXNGCZ.'1 1 , E/INCZM,  DZN*  H) 

COM  PUT  E MUTUAL  IMPEDANCE  PE.TVEEN  TU3  FA  I I. ELL  SEGMENTS 
OF  TH I " Ii.E  El  POLES  - PY  PIECE-VISE  SINUSODAL 
EX  PAX  SI  v;  I VEIGX1ING  FUNCTIONS 
COMPLEX  ZMNG,  C VP  LX,  CSCI.T 

DIME”  I-  T 'N  AC  9)  , "C9  5 , VC  9 ) , SUC9  ) , S VC  9) , Cl'C9>  , Cl  C9>  , CC9>,  S C 9>  , SPC9>, 
1SNC9),C::C9»,CPC9) 

N»  1 

CC*2. 0*C0SC  DZM) 

KF*F.*R 


C 

C 

C 
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I 


ZN"ZN  1 

DO  1 I=>1»3 

ZM=ZM 1 

DO  2 J=  1 # 3 

AC3*C  I-  1>*J)»ZN-ZM 

2 ZM*ZM+DZM 

1 ZN»ZN*DZN 

6 DO  3 1 = 1*9 

IF  CACI)  .LT.  O.O)  GO  TO  21 
UCI  ) *SQRT C fih+ AC  I ) * AC  I ) ) «•  AC  I > 

VC  I ) “HR/UC I > 

GO  TO  22 

21  VC  I >=SGhTCRR*A( I )*AC I ) )-AC  I ) 

UC  I ) =P.R/ VC  I ) 

22  CALL  SICICSI,CI.UCI>) 

SUC I ) = SI 

CUC I )=CI 

CALL  SI  Cl  C SI  # Cl  » VC  I ) ) 

SVC  I )»SI 
CVC I ) =CI 

SPC  I )«SUC  ! )«-SVC  I > 

SNC I >*SUCI )-SVC I ) 

CPC  I > = CUC I )*CVC I ) 

CNCI >=CUC I ) -CVC I > 

3 CONTINUE 

CC  1 )=*COSC AC  1 ) ) 

CC  2) »COSC AC  2) ) 

CC  3) »COSC  AC  3) ) 

CC  7) =COSC  AC  7 ) ) 

CC8)*C0SCAC8)> 

CC9)«C0SC AC9 ) ) 

SC 1)=SINCAC 1)) 

SC2)«SINCAC2)  ) 

SC  3)*SINC AC  3) ) 

SC7)»SINCAC7)) 

SC8)*SIN(AC8)) 

SC9)«SINCAC9)) 

hL»CC  1 ) »CCPC 1 ) -CPC  A ) > -SC  1 )*CSNC4)-SNC 1 ) ) ♦ C C 3 > * C CPC  3 ) -CPC  6 ) > -SC  3)* 
1 CSNC  6) -SNC  3) )*CC  7) *CCPC  7) -CPC  A) ) - SC  7 ) * C SN C 4) - SN C 7 ) ) ♦ C C 9 ) * C CPC  9 ) 
1-CPC  6) >-SC  9)*C  SNC A) -SNC 9) >-CC* CCC  2)*C  CPC  2) -CPC  5) >-SC2)*C  SNC  5> 

1 -SNC 2) )+CCR)*CCPCR)-CPC  5 ) ) - S C 8 ) * C SN  C f>)  - SN  C 8 ) ) ) 

AG=CC 1 )*CSPC 4) -SPC 1 ) >-SC 1 )*CCNC 4)-CNC 1 ) )+CC3> *C SPC 6) -SPC 3) >-SC3>* 
1 CCNC6)-C:j  f 3)  )+CC7)*CSPC4)-SPC7)>-SC7)*C  CN  C 4 ) - CNC  7 ) >+CC9)cC  SP C 6) - 
1SPC9) ) -SC  9)*  C CMC  M-CNC9) >-CC»C CC  2) *C  SPC  S)-SPC  2) ) -SC2)*CCNC  5>-CNC 
12) )*CCR)*( SPC  5) - SPC  3) )-SC8)*CCNC  5>-CNC8> ) > 

ZMNG«CMPLXChL, AG)* 1 S. / C SI NC  DZN ) *S I N C DXM) ) 

HETUhN 

END 


non 


SUBROUTINE  SI  Cl  C SI  , Cl  , X) 

Z-ABSCX) 

IF(Z-4.  >1,1,4 
l Y=(4.-Z)*(4.+Z> 

3 SI=X*( CC  C<  1. 753141E-9*Y+ 1 . 568988E- 7 > *Y+ 1 . 374 1 68 E- 5) • Y+6. 939889E4>- 
l*Y-f  1. 964882 E-2>+Y+4.  39  5509E-1  > 

CI  = ( ( 5.  772 1 56 E- 1 ALO  G( Z) ) / Z-Z* ( f < ( ( 1 . 38698 5E- 1 0*Y  + 1 . 584996E-8 > *Y 
1*1. 7257 52E-6>*Y+ 1 . 185999E-4>*Y+4.990920E-3>*Y+  1 . 3 1 5308E- 1 > > *Z 
RETURN 
T SI =SIN ( Z ) 

Y=COS( Z) 

Z=4./Z 

U*<<  C ( < ( ( (4. 048069E-3*Z-2.279 143E-2) *Z+ 5.  51 5070E-2)*Z-7. 26  1642E21- 
1*Z+4.987716E-2)*Z-3.332519E-3>*Z-2.314617E-2)#Z-1.  134958E-5)*Z 
24.6.25001  lE-2>*Z+2.  583989E- 10 
V=<  t < < c ( C < (-5. 108 699 E- 3* Z+2. 319179 E-2)*Z-6. 537283E-2>*Z 
l+7.902034E-2)*Z-4.4004  1 6E- 2) *Z- 7 . 945556E-3 > *Z+2. 6 0 1 293E-2) *Z 
2-3. 764000E-4)*Z-3. 1 224 18 E- 2> *Z-6. 64644 1 E- 7 > *Z+2. 500000E- 1 
CI=Z»< SI*V-Y*U> 

SI»-Z*(SI*U+Y*V>  +1.570796 

RETURN 

END 


SUBROUTINE  LI NVRTC Z» NV , NP, M2 > 

I MPEDANCF  MATRIX  INVERSION  ROUTINE  FOR.  PLANAR  ARR.AY 
-INPUT  SUBMATRIXES  ZO,  ?.  I » Z2»  . . . ETC. 

-RETURN  COMPACT  ADMITTANCE  SUBMAT IRXFS 
COMPLEX  ZCl) 

DIMENSION  I A(  2>  , I PC  2) 

COMPLEX  PS  1 C 50) , PS2C 50) , DFLC 1 00>  » ST(  50) » CMPLX 

I T = 07 

JT  = 08 

REWIND  IT 

REWIND  JT 


STORAGE  RE  GUI  REM ENT  ■ 


PS  1 CNP**2) 
P5P(NF*»2> 
DELC2*NP**2> 
ST(NP**2) 


CALC  DELC-1)  AND  p£(C0>,0>. 
N2=NP*NP 
DO  5 1 = 1,  N2 

DFL( I )=Z( I ) 

N=NV-  1 

CALL  LINF.r.  ( DFL.NP) 

IF(N.LF.O)  GO  TO  100 

CALL  MAT  XLT < I FL#  ZCN2+ 1),PS1»NP> 
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CALL  VI  APFC  1,N2,  PS1,  IT> 

I AC  l ) - ST AFT  ADDhFSS  OF  PSCCM-1),0) 
I AC  2 ) - STAhT  ADDKFSS  OF  FSCCM),0) 
IBC1)  - START  AEDRESS  OF  DELCM-2) 

I BC  2)  - START  ADDRESS  OF  DELCM- 1 ) 

IAC 1)=1 
I AC  2)=  1 
I BC  1 ) = 1 
IBC2)=N2+  1 
MZ=N2+N2+  1 
MM  a 0 


ITERATE  ON  M=l,2,  ...  N.  FOh  M=N,  ONLY  CALC  DELCM-l) 
DO  45  M=1,N 

IO=I AC 1)+MM 
MM-MM+N2 
I 1 ■ I AC  2) ♦MM 

10  - START  ADDRESS  OF  PSC C M- 1 ) > M- 1 ) 

11  - START  ADDhESS  OF  PSCCM),M) 


CALC  DFLCM-1) 

CALL  RTAPEC I 0»N2»PS1, iT> 

CALL  MATMLTCPS1,PS1,ST»NP) 

I J»0 

DO  20  J=  1,NP 

DO  20  I = 1 » NP 

I J*I  J*  1 

STCI J)=-STCI J) 

IFCI.EG.J)  STCIJ)*1.+STCIJ> 

CALL  LI NF.CC ST,NP) 

ID-IBC  1) 

I D 1 = I BC  2 ) 

CALL  MATMLTC  ST, DELC I D) , DELC I D 1 » , NP> 
IFCM.EQ1N)  GO  TO  SO 


CCLC  PSC  CM) , M) 

MZZ-MZ 
MS*I AC  1) 

DO  25  1 = 1, N2 

STCI )=ZCMZZ) 

MZZ-MZ L* I 
MZZ-MZ-N2 
DO  30  I S*  1 , M 

FORMAT C IX, I 20/) 

CALL  HTFFFCMS,N2»rSl,IT) 

CALL  TR.XXL  » C PS  1 , Z C MZZ ) , ST,  NP) 

ThMMLT  ACCUMULATES  IN  S7  - TR.ANSPCPSC 
MS=MS*N2 

MZZ=MZZ-N2 

MZ=MZ*N2 

CALL  MATMLTC  DELC I E I ) , ST, PS2,NP) 

CALC  PSC  C M, h)  FOR  R=0,  1,.,.  M- 1 . 


I0R»IA< 1) 

I IR-IAC2) 

I MR* I 0 

DO  40  I R*  1,M 

CALL  RTAPE<  IMR,N2»PS1»  IT) 
CALL  MATMLT (PS1»PS2,ST,NP> 
IMR*IMR-N2 

CALL  RTAPEC I 0R,N2,  PS  1 , I T) 
DO  42  1*1, N2 
PS1(I)=PS1CI)-ST(I) 
CONTINUE 

CALL  WTAPE( I 1 R,N2>  PS  1 , JT ) 

I 1R*I 1R+N2 
I OR* I 0R+N2 
CONTINUE 

CALL  VTAPECI 1,N2,PS2,  JT) 

I * JT 
JT*IT 
I T=I 
I * I A(  1) 

I A(  1 ) * I A( 2) 

I A(2)*I 
I * I B<  1 ) 

I B(  1 ) = I P(  2) 

I B( 2) * I 
M2P=NP*M2 


CALC  Y(  0, 0)  SUBMATRIX  ; Y ( 0, 0 ) = DEL C N- 1 ) 


CALL  SH I FT < Z ( 1 ) , DEL< I P< 2) ) , NP, M2) 


CALC  Yth,  0)  SUBMATRIX  ; Y( h, 0) *-PS< CN- 1 ) , ( R- 1 ) ) *DEL<N- 1 ) 

R*  1 N 


M2P=M2*NP 
I OR* I A<  l) 

I 1 R=M2P+ 1 
DO  60  1*1,  N 

CALL  RTAPE( I OR, N2#  PS  1 » IT) 

CALL  MULT 2 ( PS  1 , X ( 1),ZCI 1R),NP,M2) 
I OR* I 0K+N2 
I IK* I Ib.+  MSP 


CALC  Y(O.K)  SUPMAT  MX  J Y < 0,  r.  > * TRA.\  L ( Y ( R,  0 ) ) 
Y<  0,  fc) *- DFL<  N-  1 > * TR  AN  S ( PS  ( N - 1 >, (K-  1)  > 

R*  1, N/2 

DELCN-  1 ) *Y ( 0,  0)  = TRA.\’Sf  LFL(N-  1 ) ) 
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NM-N/2 
IOR-IA( 1) 

I lR-MPVi+l 
DO  70  I-1,NM 

CALL  RTAPEC I0R,N2, PS1, IT) 

CALL  TRANS( PS  1 » DEL( IB( 1 ) > ,NP) 

CALL  SHI  FT (ST,  DEL( IB( 1))»NP,M2) 

CALL  MULT2(DEL(IB(2)),ST,2(I  1R),NP,M2> 
I OR* I 0R+N2 
70  I 1R-I 1R+MPV 


CALC  Y(R,S)  ; S*l,...,N/2  , R=1,...,N 

Y(R,S)=Y(R- 1,S-1)*PSCR- l ) *DEL*TRANS( PS< S- 1 ) -PSCN-R)*DEL*TRANS< 
PS(N-S)) 


I J-MPW+ 1 

I 1R=(N-1)#N2>IA(  1) 

I4R=MPW* 1 
DO  90  IS*1,NM 
I J«IJ*M2P 

I2R*(N-1)<»N2+IA( 1) 

I3R-IAC 1) 

DO  80  IR«1,N 

CALL  RTAPECI  1R,N2,PS1,  IT) 

CALL  TRANS(  PS  1*  DEL( IB(  1 > >,NP> 

CALL  SHIFT(ST,DELCIBC 1>),NP,M2) 

CALL  MULT 2 ( DEL( IB(2))»ST»DEL(IB< 1)),NP,M2> 
CALL  RTAPEC I 2R,  N2»  PS  1 , I T> 

CALL  MULT2(PS1,DEL(IB< 1 ) ) , ST, NP, M2) 

CALL  hTAPEC  I3R,N2,PS1,  IT) 

CALL  MULT2(PS1,2CI4R),DEL(IBC 1>>,NP,M2> 
KK*IB< 1) 

JI*I J-(NW* 1 )*M2P 

DO  88  K= 1 , M2P 

2(IJ>*2(JI  ) ♦DEL(KK) -ST  C K ) 

JI*JI«>  1 
KK*KK*  1 
88  IJ*IJ+1 

I JT=  I J-M2P 
I3R-I3R+N2 
80  I 2R* I 2R-N2 

I IK*  I 1R-N2 
90  I 4R* I 4R+MPW 

RETURN 

100  IF(N.EC.O)  GO  TO  110 

PRINT  1010 

1010  FORM AT ( //  10X  22H I LLEGAL  CALL  TO  LINVET  //> 
110  DO  115  1*1, N2 

115  2 ( I ) = DFL< I ) 

1000  RETURN 
END 
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SUBROUTINE  MATMLT ( A.  B»  C#  NP> 

CALCULATES  C3A*B,  A,  B,  C ARE  NPXNP  MATRICES  OF 
COMPLEX  NUMBERS. 

COMPLEX  A(1)»B(1)*C(1)*D 
I J*0 
L*  1 

DO  220  Jal.NP 

DO  215  1-1, NP 

I J3I J+  1 
D30. 

KJ»L 
I K*  I 

DO  210.  K»l,NP 

D3D+AC IK)*BCKJ> 

IK3IK+NP 
210  KJ3KJ+1 

215  C( I J) *D 

220  L3L*NP 

RETURN 
END 
C 

SUBROUTINE  TRMMLT ( A»  B. C. NP) 

C ACCUMULATES  IN  C -TRANSPf A) #B 

COMPLEX  A(  1 ) . BC  1 ) # C(  1 ) # D 
I J«0 
L3  1 

DO  130  J3  1 # NP 

M»  1 

DO  125  I31#NP 
I J=I  J+  1 
Dx  0. 

KI  3M 
KJ3L 

DO  120  K=  l.NP 

D*D+Af KI )*B(KJ) 

KI-KI+l 

120  KJ=K J* 1 

M3M+NP 

125  C( I J)=C(I J)-D 
130  L=L+NP 

RETURN 
END 

C 

C GENERAL  MATRIX  INVERSION  ROUTINE 

SUBROUTINE  LINEC(C.LL) 

COMPLEX  C( 1>#ST0R,ST0,ST,S 
DIMENSION  LRC 53 ) 

COMPLEX  X 

CAPO  ( X ) =REAL  C X ) * R.F  AL(  X > ■♦■  A I MAG  ( X ) ® A I MAG  C X ) 

DO  1 A 5 I31,LL 

LRf I )*I 

1A5  CONTINUE 

M130 
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I-M.LL 


1,LL 


DO  180  M»1,LL 

K-M 

DO  155  I-M.LL 

K1«M1*I 

K2*M  1 + K 

I FC  CABQ<  C<K 1 ) ) -CABQC  CCK2) ) ) 
K»I 

CONTINUE 

LS«LR(PO 

LR(M»«LR<K> 

LROO-LS 

K2-MUK 

ST0R-CCK2) 

Jl-0 

DO  160  J-1,LL 

Kl-JUK 

K2-JUM 

ST0-CCK1) 

C(K1 )*C(K2) 

C(K2) -STO/STOR 
Jl-JULL 

CONTINUE 
K l*M  1 + M 
CCKD-I./STOR 
DO  175  I » l»LL 
I F<  I -M)  165*  175#  165 
Kl-Ml+I 
ST-CCKl) 

C(K1 >«0. 

Jl-0. 

DO  170  J»1,LL 

KI-JUI 

K2* J1*M  * 

CCK1 J-CCK1 )-C(K2)»ST 
J1-J1+LL 

CONTINUE 
CONTINUE 
Ml -MULL 

CONTINUE 

J1*0. 

DO  205  J»1,LL 
IFfJ-Lr.CJ))  185,200.  185 
LRJ»LR( J) 

J2-(LrJ- 1)*LL 

DO  195  I - 1 , LL 

K2»J2*I 
K1«J1*I 
S-C(K2) 

C( K2) =C( K 1 ) 

C<  K 1 ) *S 
CONTINUE 


155,  155,  150 


LRC  J) =LRCLRJ) 

LRCLRJ)*LRJ 

I FC  J-LRC  J)  > 183*200,  185 
200  Jl= Jl+LL 

205  CONTINUE 
RETURN 
END 

SUPROUTINE  SHI  FTC  A*  B»NP#'M2) 

COMPLEX  AC  1 >,BC  1> 

MATRIX  A(NP*M2)  IS  EXTRACTED  FROM  SQUARE  MATRIX  BCNP*NP) 
M1«=NP/M2 

K=0  _ 

L=  C M 1 / 2)  *NP 
DO  1 J=1*M2 
DO  2 I = 1 * NP 
K=K+  1 
L=L+  1 
ACK)  = BCL  > 

L“L+ ( M 1 - 1 ) #NP 

RETURN 

END 

SUPROUTINE  MULT  2 ( A,P»C*NP*M2) 

COMPLEX  AC 1 > » FC 1 > * CC 1), CMPLX*  D 
CALCULATE  MATRIX  MULTIPLICATION  OF  C=-A*P 
WHERE  A=NP<*NP  * P=«NP*M2  A C*NP*M2 
I J*0 

DO  1 J-  1 » M2 
DO  2 I = 1 * NP 
JK=  ( J-  1 ) *NP+  1 
D*CMPLX<0. 0*0.0) 

I J=I J+  1 
I K*  I 

DO  3 K*  1 * NP 
D=D*A<  I K>  * BC  JK  ) 

I K= I K+NP 
JK=  JK*  1 
CC I J) =-D 
CONTINUE 
RETURN 
END 

SUPROUTINE  TRANSC  A, D, NP) 

COMPLF.v  AC  1 ) * PC  1 > 

MATRIX  P«TRANSPJSE  OF  A 
K=  0 

DO  1 J=  1 » NP 
L*J 

DO  1 I = 1 * NP 
K=*K«-  1 
PCK)=A(L) 

L=L*NP 

RFTtlRN 

F.ND 


'-1-  vtS 


SUPi.  M’T  I*’F  RTAPFI  I 0#NP,PS,  IT) 

C JiFAP  IN  A PLOCK  OF  DATA  FROM  TAPE 

COMPLEX  PS( 50) 

REW  IND  IT 
KK*I 0-  1 
LL-KK/N2 

IF(LL.EC.O)  GO  TO  2 
DO  1 1 = 1#  LL 

1 hEADCIT) 

2 HEAD(IT)  PS 
RETURN 

END 

SUBROUTINE  WTAFFt I 0# N2# PS# I T ) 

C WRITE  ARrAY  PS( I ) ONTO  TAPE 

COMPLFX  PS<  50) 

REWIND  IT 
KK-I 0-  1 
LL=KK/NP 

IFfLL.FC.O)  GO  TO  2 
DO  1 I = 1 , LL 
1 READ(IT) 

? WRITE! IT)  PS 

hFTUKN 
END 


■ 

a 


SUBROUTINE  BLK2(M1#M2#M»A#B) 

C GENFiATE  SYMMETRIC  TOEFLITZ  MATRIX  BY  1-ST  ROW  A COLUMN 

COMPLFX  A( 1 ) , B( 1) 

Ml  2*M 1 *M2 
M123*M12*M1 
M1234=M123*M2 
DO  100  N*  1#M2 
DO  100  1 = 1 # M2 
DO  100  J«1#M1 
DO  100  L= 1 #M 1 

K 1 = 1 APSC  L- J) ♦ 1 + M 1*1 ABS( I -N  ) 

IF(I.FG.N)  GO  TO  15 
IF(I.LT.N)  GO  TO  16 

IFCL.LT.  J)  K1=K1«-M12-<I  APS!  I-N)*M1) 

GO  TO  15 

If  IF(L.GT.J)  K1=><1*M  12-!  I APSC  I -N)*M  1 ) 

1 5 K*L«-M  1 * f I - 1 ) ♦M  1 2*  C J-  1 >*M  I 23°  ! N-  1 ) + M 1 234*  ! M-  1 ) 

P! K ) = A!  K 1 ) 

100  CONTINUE 
RETURN 
END 

SUBROUTINE  NOMl'L!  M 1 # DZ#  A-C#Z2) 

C CALCULATE  IMPEDANCE  A ADMITTANCE  MATRIX  OF  A DIPOLE 

COMPLFX  Z2! 1 ) # ZMNG 
ZN-0.0 


L-  0 

DO  1 I - 1 , M 1 
ZM-0. 0 
DO  2 J=*1,M1 
L«L*  1 

Z2(L)=ZMNG{ZM,ZN,DZ,DZ, AK) 
Z2(L)=Z2CL)*0. 01 
2 ZM=ZM+  DZ 

1 ZN-ZN+  DZ 

RETURN 
END 


SUBROUTINE  FAR< SS, NR, NT, NY, NK, NS , N J, HX, HY, HZ, M, DZ, PH, TE, 
1TOT  1»T3T2,N,M1,V) 

C CALCULATE  FAR  FIELD  E-FIELD  AND  GAIN  OF  GIVEN  Z-DIRECTED 

C DIPOLES  WITH  KNOWN  CURRENT  DISTRIBUTION 

COMPLEX  AZ, TERM, YJ0»JYl»E2,SS(  1), PIN, VC  1) 

DIMENSION  PH( 1),TE( 1),T0T1( 1),T0T2(  1) 

DIMENSION  HX ( 1 ) , HY  Cl), HZ  C 1 > 

COMMON  AAC 180>,BBC 180) 

PINaCMPLXCO.,0.) 

K» ( M 1 ♦ l)/2 
DO  1 1*1, N 

PIN*PIN«-V(I  )*CONJG<SS(K)  ) 

1 K«K*M1 
NA*NR+ 1 

PO WFE=REAL (PIN) 

NB-NT+ 1 
NC=NK*  1 
ND-NS*  1 
LL*0 
PIG-0. 0 

DO  11  LX-NC,  ND,  N J 
LM«LX-  1* 

THE- 0. 0 1 7 A 53 3* FLO AT (LX- 1 ) 

CTH-COSC  THE) 

STH-SI N C THE) 

STST=STH*STH 

K2*0 

DO  11  LK*NA»NB,NY 
LO-LK- 1 

PC»0. 0 174533*FL0ATCLK- 1 ) 

PC-COSf  PC) 

PS-S INf PC) 

AZ*CM?LXf 0. 0, 0. 0) 

DO  4'  A J=  1 , N 

YZ*  STH® ( PC’HX  <J)+HY(J)*PS) 

YJ0=CMPLX<  0. 0,YZ) 

TERM  = CF.XP(  YJO) 

TS-SINCDZ) 

TC-COSC  DZ) 
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MOC  J-  1)*M1»  1 
MB* J*M 1 
ZDZ*HZ< J) 

DO  47  I *MC#MB 
ZDZ*ZDZ+DZ 

YJ1*CMPLX( 0. 0, ZDZ#CTH) 

E2*CEXP( YJ1 ) 

IFILX.EQ.  1.0R.LX.EQ.  181)  GO  TO  46 

AZ*AZ*2.*TERM*E2/<TS*STST)*( C0S< DZ*CTH> -TC)*SS( I > 
GO  TO  47 

46  AZ*AZ+TERM*E2*DZ*SS(I) 

47  CONTINUE 
44  CONTINUE 

LL*LL«- 1 

TT«CABS(AZ)*STH 
IF(TT.GT.BIG)  BIG*TT 

PP*ATAN2< AIMAG<  AZ)#REAL<  AZ) )/0»  0174533 
PP*(  377./ 12«  5663706) *TT** 2/ POWER 
AA(LL) *TT 
BBCLLJ-PP 

TE(LL) *TE(LL)+TT**2 
TOT1 (LL) “TOT  1 (LL)*TT 
PH ( LL  > * PH ( LL ) ♦ PP*  * 2 
T0T2ILL) *T0T2( LL)+PP 
11  CONTINUE 
M-LL 

92  RETURN 
END 


SUBROUTINE  STDV( TOT  1 # TOT  2# PH#TE#Nh»NT» NY»NK#  NS#N  J# ST1#  ST2# IN) 
C DATA  PROCESSING  OF  THE  FAh  FIELD  POINTS 

DIMENSION  ST  1 ( 1 ) # ST2( 1 ) # TOT  1 < 1 ) # T0T2( 1)#  TE< 1 ) # PHC I ) 

PRINT  13 

13  FORMAT C///8X,  'THE'.  7X#  'PHI  '#6X,  'MEAN  E-FIELD'# 

SISX.'STD  '#  1 OX#  'MEAN  GAIN '#  16X#  'STD'#///) 

NA-NR* 1 
NB-NTM 
NC*NK» 1 
ND*NS-»  1 
LL-0 
INI-IN-1 

IFCIN1.EQ.0)  INI-1 

DO  11  LX*NC»ND#NJ 
DO  11  LK=NA#NB#NY 
LL-LL*  1 

ST1<LL)=*SQP.TC  APS( TEfLL ) - TOT  If  LL ) » * 2/ FLOA T C IN)  ) /FLOAT*  INI)) 
ST2(LL>  =bGHT* APS* PH< LL > -T0T2C LL) **2/rL0 AT ( IN) )/FLQAT< IN  1 ) ) 
TOTl(LL)=TOTl( LL)/ FLOAT ( IN) 

T0T2<LL)*T0T2(LL)/rL0AT< IN) 

ITH-LX-1 
IPC*LK-  l 
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PRINT  12. ITH, IPC,T0T1CLL>,  ST  1 C LL) , T0T2C LL) , ST2C LL) 
11  CONTINUE 
MAX=CLL* 1 )/2 

VRITEC2,  15)  T0T2C  MAX) , ST2CMAX) 

PRINT  15*  T0T2C MAX ) , ST2C MAX ) 

> FORMATC/// IX,  ' AVERAGE  BROADSIDE  GAIN  = ',F20.6/, 
*1X,  ' STD  = ',F20.6/) 

» FORMATC 1X,2I 10, 1P4E20. 6) 

RETURN 

END 


SUBROUTINE  TAPERC  V, M 1 » M2, M3, HX, HY , HZ) 

NON-UNIFORM  EXCITATION  MULTIPLIER 
COMPLEX  V( 1) 

READ  50#  RX,  RZ 
.0  FORMATC  5F 10.  5) 

DIMENSION  VVC iOO),PPC 100),HXC 1),HYC 1>,HZC 1) 

PRINT  77 

7 FORMATC// IX#  'MAGNITUDE  TAPER  * PHASE  SHIFT  IN  EACH  DIPOLE'#//) 
M23*M2*M3 

CALL  TAPCM2,M3,HX,HZ,  VV) 

CALL  STEERC M2, M3#  RX# RZ#  HX#  HZ,  PP) 

PRINT  40 

lO  FORMATC// IX,  ' ELEMENT  '#  5X#  'MAGN I TUDE  TAPER  4 , 5X, 

5, 'FHASE  SHIFT',//) 

PRINT  10#  CI,VVCI)#PPCI  )#  1 = 1, M23) 

DO  1 1 = 1 ,*M23 

0 FORMATC IX, I 10, 2F17. 5) 

VC  I )=VC I )*VVC I )*CMPLXCCOSCPPCI ) ), SINCPPC I ) ) ) 

1 CONTINUE' 

RETURN 

END 

SUBROUTINE  STEERC M2#  M3#  RX,  RZ,  HX,  HZ,  PP) 

PROGRESSIVE  PHASE  SHIFT  ON  EACH  DIPOLE  - 
STEERING  THE  MAIN  BEAM  IN  POTH  DIRECTION 
DIMENSION  PPC 1)»HXC 1),HZC 1) 

PI2=6. 2831853 

THX*  C 9 0, -RX)* PI  2/360# 

THZ=C90.-RZ)*PI2/360. 

L=  0 

XK=  C HX C ? ) -HX C 1 ) )*COSCTHX) 

ZK=  C HZ  C X?>  1 ) -HZC  1 ) ) * CO  SC THZ ) 

DO  100  1=1, M3 
DO  100  J= 1 , M2 

L=L*  1 

PPC  L ) = FLO  Al  C J-  1 ) * XX*  FLO  AT  C I - 1 ) * Z K 
100  CONTINUE 
RETURN 
END 


mmmm 


SUBROUTINE  TAPCM2, M3,  HX,  HZ,  VV> 

MAGNITUDE  TAPER  OF  EXCITATION  IN  BOTH  DIRECTION 
DIMENSION  HX(  1 > , HZ  C 1 ) , W<  I) 

PI2=6. 2831853 

HFX=  ( HXC  M2> -HX ( 1)>*0.5 

HFZ=(HZ(M2*CM3- 1)+ 1>-HZ< 1))*0.5 

L=0 

DO  100  1 = 1, M3 
DO  100  J=  1 , M2 
L*L>  1 

X=  C HX(  L > -HFX ) /PI  2 
Z=(HZCL)-HFZ)/PI2 
VV( L ) = FUNX ( X > * FUNZ ( Z ) 

CONTINUE 

hETUhN 

END 

FUNCTION  FUNX(X> 

FUNX»FXP< - ( ABS(X) ) > 

RETURN 

END 

FUNCTION  FUNZ(Z) 

FUNZ  = F.XP<-<  ABSCZ)  ) ) 

RETURN 

END 


SUBROUTINE  MULT ( Z, V, CUR, M 1 , M2,  M3 ) 

GIVFN  HALF  OF  THE  COMPACT  Y-MATRIX  * VOLTAGE,  CALCULATE  CUR. KEN 7 

COMPLFX  Z< 1),VC !),CURC 1),A 

M123=M1*M2*M3 

M23=M2»M3 

N2=(M23*.l  )/2 

IF(M3.LT.M?)  CALL  OKTFhf V, 1 , M2, M3) 

DO  1 I = 1 , M 1 23 
A*CKPLX( 0. 0, 0. 0) 

I 2 = M 1 23- I ♦ 1 

DO  2 J*  1 »N2 

J2=M23- J* 1 

L 1 * ( J- 1 )*M123M 

L2»C  J-  ! )<*M123+I2 

IF(J.NF.J2)  A«A+Z(L2)*V< J?) 

A*A*Z<L1 )*V< J) 

CONTINUE 
CUr.f  I ) * A 
CONTINUE 

IF(M3.LT.M2>  GO  TO  3 
RETURN 

CALL  OrlFRCClL  , Ml, M3, M2) 

CALL  ORDER <V, 1,M3,M2) 

RETURN 

FND 
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SUBROUTINE  MULT3(  Z2»  V,  CUR2, M 1 , M2, M3) 

COMPLFX  V( 1>,CUR2< 1>,Z2< 1) 

CALCULATE  THE  CURE  ENT  ON  EACH  F.LEMFNT  OF  AN  ARRAY 

WITHOUT  TAKING  ACCOUNT  THE  MUTUAL  EFFFCT 

M23=M2*M3 

MM=< (M1+  l)/2-  1)*M1 

L=  0 

DO  1 I 3 1 * M23 
DO  2 J=1,M1 
L=LM 

CUR2(  L ) =Z 2f  MM-f  J)  *V(  I ) 

CONTINUE 

CONTINUE 

RETURN 

END 


SUBROUTINE  0 RDER< V, M 1 , M2, M3 > 

REORDER  THE  SEQUENCE  OF  EACH  DIPOLE  IN  ARRAY 
COMPLEX  VC  1),VV(300) 

L=  0 

DO  1 1 = 1,  M2 
K*  I 

DO  2 J= 1 i M3 
N=CK- 1 >»M1 
DO  A M=  1 , M 1 
N=N+  1 
L=L*  1 

VVCL)=VCN) 

CONTINUE 
K=K*M2 
CONTINUE 
M123=M1*M2*M3 
DO  3 1 = 1,  M 123 
V( I)=VVC I ) 

RETURN 

END 


susnc  utine  vcltcv,vp,m2,:i3,  iy*Ai-n,AH2*  sdi,sd2> 
gz:j  spate  random  eppop  excitations  fop  appay- 

NOP  ilAL  DISTRIBUTION  OF  EOTH  MAGNITUDE  £ PHASE 
EPP.  CP  APE  ASSUMED 
COM  FLEX  VC  1 >,VPC  1 ) 

M23=M2*M3 

DO  200  1=1, M23 

CALL  PUDCGl,Atll,  SD1,  IX) 

I F ( G 1 ) 2,  3,  3 

CALL  P.UDC  G2,  AI12,  5D2,  IX) 

I F ( G 2-3*14159265) A, 5, 3 
IFCG  2+3*14  159265)3,5,5 
VC  I )=  w(I ) *G  1 *CMPLX  (C0SCG2),  S1MCG2)  ) 

CO.'J  TIMUE 
r ETUPH 
END 


D,  IX) 


SUED  OUT  I ME  F.NDCG,  AMEAM, 
A=  0 * 0 


I F ( S"D.  GT  * 3 • 14  1593)  KD=  1 

DO  7 I = 1 , KD 

IY=IX*65539 

IFCIY)  10,12,12 

I Y=  - IY 

Y FL= IY 

IX  = IY 

Y FL=YFL*  * 29  I C3C3C5E-  1 0 
A=A+YFL 

G=<A-6*0)*STD+AMEAN 
I FCSTD.GT.3* 141593)  G=<A-0 
P.ETUPN 


